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Introduction 



In this paper we introduce r-functions for the Dirac operator on the Poincare disk 
based on the formahsm in [12]. The formahsm of the first author in [11,12] is in 
turn a geometric reworking of the analysis of Sato, Miwa, and Jimbo [15] (SMJ 
henceforth) . The prototypical examples of the holonomic fields that are the central 
objects of the SMJ theory are the scaling limits of the two dimensional Ising model 
from above and below the critical point. The SMJ theory of holonomic fields pro- 
vides a beautiful setting for the earlier result of Wu, McCoy, Tracy, and Barouch 
[1,8,17,20] that the scaled two point functions of the two dimensional Ising model 
can be expressed in terms of Painleve functions of the third kind. Our work also 
owes something to the papers of [18,6] on monodromy preserving deformation the- 
ory in the Poincare disk. In the course of explaining our work we will point out 
differences with [18]. 

Our principal goal in this paper is to define and analyze the correlation func- 
tions for a family of quantum field theories on the Poincare disk that are analogues 
of the holonomic quantum field theories defined in the Euclidean plane. Below we 
will sketch how the theory goes in the Euclidean version in [12] and we will point 
the reader to the parallel developments in the hyperbolic case. 

Our starting point is a formula for the Schwinger functions of holonomic fields 
that can be written 

T(ai,a2,...,a„) = det(D'^'^), (1) 

where T(ai, a2, . . . , a^) is the vacuum expectation for a product of quantum fields 
$A (flj) [11]- In this formula the operator D"-'^ is a singular Dirac operator (on Eu- 
clidean space) with a domain that incorporates functions with prescribed branching 
at the points aj. The functions in the domain of L)"''^ are to have monodromy mul- 
tiplier e*^'^'^^ in a counterclockwise circuit of aj. In [12] formula (1) is given rigorous 
mathematical sense by first localizing D"^'^ away from branch cuts emerging from 
the points aj. The localization for D"''^ takes place in the exterior of a set S which 
is the union of strips containing the branch cuts. The localization of W^'^ is com- 
pletely characterized by a family of subspaces W"'''^ C {dS) of the Sobolev space 
on the boundary, dS, of S. These subspaces are the restriction to dS of functions 
which are "locally" in the null space of D"^'^. The subspaces W°''^ are contained 
in a restricted Grassmannian of subspaces of H'2{dS) over which there sits a de- 
terminant line bundle [16]. As is explained in more detail in [11,12], defining a 
determinant for the family D"'-^ is morally equivalent to finding a trivialization 5 
for the det* bundle over the family of subspaces W"'''^. Comparing the canonical 
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section a of det* with 5 gives us a formula for the determinant, 

The analogue of this formula on the Poincare disk is (6.7) below. Both the canonical 
section, a, and the trivialization, 5, are determined by (different) projections on the 
subspaces W°-'^. An explicit representation of these projections is obtained in terms 
of the Green function G""'^ for D"-'^ . The analogue of this projection formula on 
the Poincare disk is (6.2) below. In order to analyze (2) it is useful to compute the 
logarithmic derivative 

da log r. 

The resulting formula simplifies dramatically because the derivative of the appro- 
priate projection in the a variables is finite rank. This in turn is a consequence of 
the analogous result for the Green function, G'*''^, namely. 



da^G^'\x,y) = c,WAx)^W:{y). 

For simplicity we have taken some liberties in displaying only a schematic version 
of (3). The precise version of (3) in the Poincare disk is (4.50) below. The func- 
tions W^{x) and W*{x) are multivalued solutions to the Dirac equation which are 
branched at the points a^, j = 1, 2, . . . , n with monodromy multipliers e^'^*'^^ . Be- 
cause of their role in (3) we refer to the wave functions Wi,{x) and W*{x) as response 
functions. The formula for d log r which one obtains is 

n 

d log T = ^ {Ajyda^ + B^ddu] , (4) 

where Ai, and B^, are low order Fourier coefficients of the response functions Wi, 
and W* in the local expansions of these functions about the points a^. The version 
of this formula in the Poincare disk is (6.8) below. 

We assemble the 2n response functions into a single column vector 



W(a;) := 

with 



W{x) 
W*{x) 



W{x) = [W^{x),W2{x),---,Wn{x)\ 
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and 




d^,a'W{x) = n{x, a)W{x). 



(5) 



in the variables {x,a), where Q is a matrix valued one-form whose entries are deter- 
mined by the low order local expansion coefficients for Wi, and W* near the branch 
points ttj. The analogue of (5) in the Poincare disk is contained in (5.3) and (5.20) 
below. The consistency condition for the holonomic system (5) is the zero curvature 
condition 



These are called the deformation equations. The analogue in the Poincare disk is 
(5.26) below. These deformation equations characterize certain of the low order 
expansion coefficients of W^^ and W*. In particular it is possible to express the 
coefficients A^, and B^, that occur in the formula for the logarithmic derivative of 
the T— function in terms of the solution to (6). The version of this relationship in 
hyperbolic space can be found in (6.16) and (6.17) below. The existence of the 
holonomic system and the resultant deformation equations is, incidently, the one 
part of this theory that depends strongly on the fact that the Poincare disk is a 
symmetric space. The global symmetries of Dirac operator on the Poincare disk 
give rise to the holonomic system for W. 

In conclusion one sees that it is possible to characterize the logarithmic deriva- 
tive of the T— function in terms of certain solutions to the deformation equations. In 
the special case of the two point function in the Euclidean problem SMJ showed that 
it is possible to integrate the deformation equations explicitly in terms of Painleve 
transcendents of the fifth kind. In the special Ising case a reduction to Painleve 
III is possible (see [8]). In §5 and §6 of this paper we carry out the analysis of the 
two point case in the hyperbolic setting. We find that the two point deformation 
equations can be integrated in terms of Painleve transcendents of type VI (Theorem 
5.0) and that the logarithmic derivative of the r— function is a rational function of 
the solution to Pyi and its derivative (see (6.19)). 

We can now explain the principal differences between this paper and [18]. First, 
because the response functions are fundamental in this approach to holonomic fields, 
the basic finite dimensional space of wave functions for us is the 2n dimensional 
space spanned by W^, and W*. In [18], as in SMJ [15], the basic finite dimensional 
space is the n dimensional space spanned by the wave functions that are globally in 
L^. Formulating the holonomic system in this setting leads to extra complications, 
and the deformation equations are not explicitly worked out in [18]. Because of 
the difficulties with a direct formulation of the deformation equations, one of us 



c/0 = O A O. 



(6) 
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introduced a hyperbolic Laplace transform of the associated holonomic system and 
found that the result was a Fuchsian system of ordinary differential equations in 
the complex plane. Tau functions for this transformed system were introduced by 
adopting the r— functions for the associated Schlesinger theory (see also [15]). In [6] 
the behavior of these r— functions was examined in the limit R ^ oo and shown to 
coincide with the r-functions for the Euclidean theory. We have yet to explore the 
relation between the r— functions introduced in [6,18] and the ones we consider in 
this paper. The algebraic complexity of both sets of deformation equations makes 
comparisons difficult. 

This last point suggests an important question. Why should one believe that 
the T— functions we have introduced in this paper are the Schwinger functions for 
some quantum field theory on the Poincare disk? The reason we believe that there 
is such a relation is to be found in [12]. In that paper a transfer matrix formalism 
was established that allowed one to connect formula (2) above with the scaling limit 
of some lattice field theories defined in [10]. It is also this connection that suggested 
the "minimal singularity" definition for the domain of D°^'^ that we adopt here. We 
do not think that it would be economical to try to produce a lattice theory which 
would scale to the theory we present here. Because of the connection with the Ising 
model, however, this is an interesting problem. We should remark at this point that 
the analogue of the Ising model on the Poincare disk is not considered in this paper. 
There is such an analogue and the deformation theory in §5 is relevant to its analysis, 
but the singular Dirac operator whose determinant gives the correlation functions 
(for the T < Tc scaling limit) is different from any of the minimal singularity theories 
for Xj — ±1/2. We expect to discuss this analogue in another place. To return to 
the quantum field theory connection; there is a preprint of Segal on conformal field 
theory which sketches a "vertex operator" formalism that appears well suited to 
address such questions as Osterwalder-Schrader positivity. One of us is currently 
investigating this formalism. 

Finally we would like to mention that this paper is organized in almost precisely 
the reverse of the logical progression (1) — (6) above. The reason for this is that 
the existence theory for the response functions and the Green function proceeds 
most simply in such reverse order. We have included a table of contents to help the 
reader obtain some perspective. 

The authors would like to express their gratitude to R. Narayanan whose notes 
on local expansions on the disk were very helpful in the early stages of this work. The 
second two authors wish to thank T. Miwa for his hospitality and support at RIMS, 
Kyoto University and to the National Science Foundation whose support through 
grants DMS-9001794, INT-9106953, and DMS-9303413 is gratefully acknowledged. 
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1 The Dirac Operator on the Hyperbohc Disk 



A covering of the frame bundle. We will begin by identifying a Dirac operator 
on the hyperbolic disk, T)r, of radius R. The Dirac operator will be defined in the 
usual fashion except that we will use the simply connected covering of S0(2) rather 
than the two fold covering in our construction. This is done principally because 
the extra parameter k associated with the representation theory of the reals does 
not significantly complicate any of our analysis and we believe its inclusion might 
prove illuminating. Incidently, the radius R is carried along in all the subsequent 
calculations precisely because the zero curvature limit J? ^ oo is interesting [6]. 

The additive group of real numbers, R, is the simply connected covering group 
of the special orthogonal group, SO (2), in two dimensions with the familiar covering 
map 

cos A — sin A 
sin A cos A 



R 9 A ^ R{\) 



e S0(2). 



To emphasize the role R plays as a covering group for SO (2) we will sometimes 
write 

R= S0(2). 

Next we define a principal SO (2) bundle on Djj, which we call S(Dji) together with 
a bundle map, 

TT : ^(Dh) ^ F(Dfi), 

onto the frame bundle, F{T)f>), that is compatible with the map R. If (7 G S0(2) 
and s sg is the right action of SO (2) on the bundle S(Dr) then the compatibility 
we refer to is, 

7r{sg) = n{s)R{g), (1.1) 

where R{g) G SO (2) acts on the frame bundle in the usual fashion. By fixing a 
trivialization of the frame bundle for Y)r we will be able to use (1.1) to define the 
bundle S(Dij). The metric we wish to consider on Dji is, 

{R-^-xl-xlf 
with the associated invariant measure 

R'^dxidx2 

dit := o- 

(i?2 - a;? - a;i)^ 
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We fix an orthonormal frame, 



ei := 1 



xl + d 



dxi 



62 := 1 



™2 1^ ™2 
1 * 2 

i?2 
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which identifies the frame bundle, F(Dij) with 



r>R X so(2). 

Relative to this trivialization of the frame bundle we identify S^Dn) = D/j x S0(2) 
and the map n with the covering given by 

X S0(2) 3 {p, A) ^ 7r{p, A) = (p, i?(A)) e x S0(2) 



for p G Dij. There is a covariance of this construction with respect to global rigid 
motions of D/j that it will be useful for us to understand. Since the (orientation 
preserving) rigid motions of Djj are simplest to understand as fractional linear 
transformations we will introduce the complex notation z = xi + 1x2- Consider the 
fractional linear transformation 



w 
R 



(1.4) 



This will be an orientation preserving rigid motion of D^^ provided that 



A B 
B A 



eSU(l,l), 



where g e SU(1, 1) if and only if 



"1 


" 




"1 " 





-1 


9 = 


-1 



and det{g) = 1. The action of a transformation (1.4) on the frame [ei, 62] is simplest 
in the complexified tangent space and so we also introduce the complex vector fields 

e := i (ei - ie2) 
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and 



For g = 



A B 
B A 



write 



e := I (ei + 162) ■ 



^ ^ A^ + B 

giz) = 

+ ^ 



For g G SU(1, 1) the pair (gr, dg) induces a transformation of the frame bundle which 
we calculate in the trivialization {61,62} by first determining the action of dg on 
the complex vector fields e and e. One finds 



dgz{e{z)) = Ug{z)e{g{z)), 



and 
where 



dgz{e{z)) = Ug{z)e{g{z)), 



u„(z) := =^ — = 
^ B^ + A 



(1.5) 
(1.6) 
(1.7) 



is a number of absolute value 1. Now write 



Ug{z) + U g(z) iiUg{z) - Ug{z) ) 

-iiUg{z) - Ug{z)) Ug{z) + Ug{z) 

Then without difficulty one sees that (1.5) and (1.6) imply 

{dggei{z),dgge2iz)) = {ei(g(z),e2igiz)))g'(z). 



(1.8) 



So that in the trivialization of the frame bundle determined by {61,62} the map 
induced by {g, dg) is 

iz,r)^ig{z),g'iz)r), (1.9) 

where r e SO (2). Next we will show that the map (1.9) is covered by a map on 
SCDr). For fixed g e SU(1, 1) this clearly comes down to defining a continuous 
logarithm for D/j 3 z ^ U'gi^) ■ Since the map Ug{z) is continuous (|A|^ — = 1 
for g G SU(1, 1) implies that Bj^ + A is never for z G T)r) it follows that the 
image of the map Ug{z) is simply connected and one can find a continuous map 
z ^g{z) SO that 



Ug{z) 



The bundle map. 



{z,X)^{giz),X + Xg{z)), 
7 



(1.10) 
(1.11) 



on X SO (2) then clearly covers the the transformation (1.9) on the frame bundle. 
Of course, the map Xg{z) is not unique; for any integer n the map Xg{z) + 2n7r would 
serve as well. The fundamental group of SU(1,1) is Z and so we do not expect to 
make a canonical choice of such a logarithm for all elements of SU(1,1). We do 
expect that the group of bundle transformations of S{Y)fi) which covers the action of 
SU(1,1) on F{Y)fi) is the simply connected covering group SU(1, 1). However, since 
there is no simple (i.e. matrix) model for SU(1, 1) and since we are more interested 
in the lifts of individual elements and of one parameter subgroups in SU(1,1) we 
will not pursue this matter. However, it will be convenient for descriptive purposes 
to refer to a pair {g, Xg{z)) where satisfies (1.10) as an element of SU(1, 1). 

The Dirac operator. We will now introduce a Dirac operator and observe that 
the lift of the action of an element in SU(1,1) to the bundle SCDr) also acts on the 
solution space of the Dirac equation. 

The Dirac operator we define will act in the space of sections of an associated 
bundle for 8(0 r). The bundle we have in mind is one associated with a complex 
unitary representation of SO (2) given by 

SO(2)9 A^rfc(A) := 

where k is an arbitrary real number. We write 







■<fe-|)A 



(1.12) 



for this associated bundle. One can view the elements of this bundle in the usual 
fashion as equivalence classes of pairs (s, v) where s e SCDr) and v e and (s, v) 
is equivalent to {sX,rk{X)v) for any A G S0(2). Because we are working with a 
fixed trivialization of the bundle SCDr), the associated bundle SCDr) Xfc can 
be trivialized by the map 

{z,X,v) ^ {z,rk{X)-\) e Dfl X C^. 

To avoid confusion we will always work in this trivialization of the associated bundle. 

The reason for choosing the representation (1.12) to define a bundle of spinors 
for the Dirac operator has to do with the existence of a complex representation of the 
Clifford algebra of on which behaves well with respect to the representation 
(1.12). The representation is generated by 



71 



1 

1 
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and 



If we define 



72 



-i 

1 



7(^) := 71^1 + 72^^2, 
then one can easily confirm that 



(1.13) 



This is the relation that allows one to ^fine a Dirac operator that will have the 
appropriate covariance with respect to SU(1, 1). 

In order to write down the Dirac operator we first recall some facts about 
connections. On a principal bundle, once a local frame has been picked, a connection 
is represented by a Lie algebra valued one form over the domain of the trivialization 
in the base. In our case we start with the global frame {61,62}. The Levi-Civita 
connection form u for this frame (which is determined by the torsion condition 

Vei62 - Ve2 6i - [ci, 62] = 0) is 



and 



•^(62) 



2X2 
i?2 



2x1 



-1 

1 



-1 

1 



Since the map it : SCDji) — > F(Dij) is a covering map on the fiber one can pull back 
the connection on F(Djj) to a connection, Q, on 5'(Djj). One finds the connection 
one form relative to our standard trivialization of SCDn) is 



a;(6i) = 



C^(62) 



_2X2 



where the Lie algebra of SO (2) is identified with R in the usual fashion. The 
connection cD on the principal bundle 8(1^^) induces a connection on the associated 
bundle S{T)r) X/. C^. In the trivialization of the associated bundle defined above 
one finds the connection one form 



drku){ei) 











2X2 
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drku){e2) = 



i{k+\) 
i{k-\) 



2x1 



The Dirac operator (in the triviahzation Tin x defined above for the associated 
bundle) is then 



Dk = XI 7j {ej + drkui{ej)) 



(1.14) 



or in complex notation 



Dk = 





k- 



where 



and 



fe-T 



:= 2 1 



i?2 



{k 



-2 



1 - 



i?2 



(fc+l) 

i?2 



(1.15) 



(1.16) 



We have written for the formal adjoint of with respect to the inner product, 



f{x)g{x)dn{x), 



on D/j X C. Thus the Dirac operator (1.14) is formally skew symmetric with respect 
to the inner product, 

(1.17) 



/ 



f{x) ■ g{x)dii{x), 



on D^^ X where u ■ v = uiVi + U2V2- The indexing by k in (1.14) and (1.15) has 
been done so that /c = corresponds to the Dirac operator that one would obtain 
from the usual two fold covering of SO (2) by itself given by 

9^9'. 

We remark that i^T^ is the familar Maass operator (see, e.g., [2]). 

Covariance of the Dirac operator. One of the principal objects of study for 
us will be multivalued solutions to the Dirac equation -DjtV' = (or more precisely, 
solutions to the massive version of the Dirac equation {ml — Dk) •0 = 0) that are 
branched at points a-,- e Hn for j = 1, 2, . . . n. In order to understand these solutions 
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to the Dirac equation in a neighborhood of the point aj it wiU be convenient to 
employ coordinates for Tin which are centered at aj. Of course, we can use a 
fractional linear map, 



w = R- 



Z — Gn 



ajz 



'i/j{z{w)) no 



to make such a change of coordinates. However, the function w 
longer satisfies the Dirac equation in the w variables. This inconvenience^an be 
overcome by using the covariance of solutions to the Dirac equation under SU(1, 1) 
that is a consequence of the fact that the Dirac operator is a natural operator on 
the associated bundle of spinors which depends on a choice of spin structure and on 
the representation of the spin group but not on the choice of local frame or base 
coordinates used to realize (1.14) as a differential operator. We will now explain 
this covariance from the passive point of view in which fractional linear maps on 
IDji are regarded as a change of coordinates on Djj. Let 



w 



Az + BR 

' = qiz) = R= 

^ Bz + AR 

denote a fractional linear map with 



A B 
B A 



e SU(1,1). 



The inverse coordinate change is 



z{w) = R^ 



Aw + BR 



Bw-AR 



Also 



where 



eiz) := 1 



i?2 



dz = v{w) 1 



\w\ 



v{w) 



i?2 

Bw - AR 
Bw-AR' 



Taking complex conjugates one has 



dyj = v{w)e{w) 



e{z) = v{w)e{w) 

as well. Thus the coordinate representation of the Dirac operator in the w coordi- 
nates is 





v{w)e{w) + ^-^T^z{w) 



v{w)e{w) - '"^j^i\ {w) 
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Now choose a continuous logarithm, iX{w), for the function w — > v{w) so that 
and define 



:= e^'\ 



for I e R. Let 



yl^ + h 
f''"^ 



Then the covariance property of the Dirac operator under the coordinate change 
z ^ w is simply expressed by 





e[w) -\ -^^r^w 



e{w) 



w 







(1.18) 



Observe that the right hand side of (1-18) has precisely the same form in the w 
coordinates that the Dirac operator has in the z coordinates. Thus if iIj{z) is a 
solution to the Dirac equation Djfip = in the z coordinates, then M^^(w)ip{z(w)) 
will be a solution to exactly the same equation in the w coordinates. An observation 
which is useful in checking (1.18) with a minimum of computation is 



with 



e{w) (v) 



e{w) (v'') = lv^e{w) (iX) = lv''v ^e{w) (v) , 



i?2 



Bw-ARj 



B 



B? J \Bw - AR 



and 



with 



e{w) (y) = Iv e(w) (iX) = —Iv ve(w) (y ^) , 



/ \ / i\ / . I^T \ I Bw — AR\ \w\ 
e{w) (v-^) =1-^1 a,„ I = — 1 = 11 



i22 



Bw-AR 



B 



R^ J \ Bw-AR J ' 



We can rephrase this observation in a slightly different form that will be useful for 
us: 
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Proposition 1.0. Suppose that \E'(2;) is a solution to the Dirac equation 

(m-Dfe)* = 0. 



If 9-- 
that 



A B 
B A 



e SU(1,1) and z ^g{z) is a continuous function of z defined so 



then the function, 

z ^{z) = 



v{g,z)'^'^^ 

v{g,z)'^~2 



remains a solution to the same equation (m — -Dfc)* = 0. 



(1.20) 



It will also be useful to identify the infinitesimal symmetries of theJDirac operator 
associated with lifts of the one parameter subgroups of SU(1,1) into SU(1, 1). It will 
be convenient to introduce the notation, 

cht = cosh t, 
sht = sinh t, 
tht = tanh t . 

One choice of three one parameter subgroups of SU(1,1) that generate the Lie 
algebra su(l,l) is given by 

ch t sh t 
sht cht 



t 



with generator 



and 



with generator, 



and 



1 

1 



ch t ish t 
—ish t cht 



X2 = 



i 
-i 



e** 
e-** 
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with generator, 
For 

one finds 



X, 



i 
-i 



9j{t) ■■= e 



-tXi 



v{g2it),z) = 



and 



R-th {t)z 
R-th {t)z '■ 

R - ith {t)z 
R + ith {t)z 

2it 



v{g^{t),z) = e'^\ 



For j = 1, 2 we choose the function \g.(^£){z) by normahzing it so that \g.(t){Qi) = 0. 
For J = 3 we choose )^g^{t){z ) = 2t. Define 



F,* := — 

^ dt 



v{gj{t),z) 






v{gj{t),zy 



(1.21) 



t=o 



Then the Fj are infinitesimal symmetries for the Dirac operator. It is very 
convenient to introduce the complexified infinitesimal symmetries, 



Ml := i (Fi - iF2) , 
M2 := i(Fi + zF2), 



and 



(1.22) 
(1.23) 

(1.24) 



The reason this is convenient is that the commutation relations for Mj are 



[Mi,M2] = -2M3, 
[Mi,M3]=Mi, 
[M2,Ms] = -M2. 

The last two show that Mi and M2 are ladder operators on the eigenspaces of 
M3 and this will allow us to parametrize local expansions of solutions to the Dirac 
equation in a sensible fashion. 
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§2 Local Expansions 



Eigenfunctions for infinitesimal rotations. We next wish to consider multival- 
ued solutions to the Dirac equation, 

which are branched at a point a e D/j with fixed monodromy about a. Following 
SMJ [15] our principal tool in understanding such solutions will be to look at ex- 
pansions in eigenfunctions for the infinitesimal "k-rotations" about the point a. To 
begin we consider the case a = 0. At a = the infinitesimal generator of k-rotations 
is M3. In geodesic polar coordinates (r, 6) centered at we have 



-ids + + I 



-ide + /c - I 



where the relation between the complex Cartesian coordinate z = xi + 1x2 and 
geodesic polar coordinates is 

z = Re'hh ^ . (2.0) 

Of course, geodesic polar coordinates are not uniquely determined at a = 0. One 
must also specify the geodesic ^ = 0. Thus when we speak of geodesic polar coor- 
dinates we will suppose that a choice of the geodesic 9 — has been made. This 
will be important for us since the local constructions in this section will give rise 
to multivalued functions. These multivalued functions will have different branches 
and in this section we will pick out such branches by using the geodesic ray ^ = 
as a branch cut. 

We turn now to the eigenvalue equation, 

M3^{r,9)^{k + X)^{r,e), (2.1) 

whose solutions will have a monodromy multiplier, 

after a counterclockwise circuit of 0. We will now analyse the multivalued solutions 
of the Dirac equation (m — D^)^ = which are also eigenfunctions (2.1) of M3. In 
studying solutions to the Dirac equation (m — D^)^ = it is very convenient to 
introduce the associated Helmholtz equation 

(m- Dk)*im- Dk)'i? = (rri^ - Dl)^ = 0. (2.2) 
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In the notation of the previous section this is 



We introduce 



A, = -KK, + 
and a httle calculation shows 

A«+i = -K^K* + 



Ak{k + 1) 



* = 0. 



4k{k + 1) 

^2 



For reasons that will be clear later on we introduce a parameter s by 



(2.3) 



(2.4) 



(2.5) 



There are two roots to (2.5) and for definiteness we will in the future restrict our 
attention to the solution s of (2.5) with Rs > 1. Another form of (2.5) which is 
useful for us is 

2\ 2 4/fc2-i 



If one rewrites the Helmholtz equation using (2.3)-(2.5) the result is 







* = s I s - - 1 



(2.6) 



To analyze the solutions of the Dirac equation which are eigenfunctions of M3 above 
we will first consider the problem of finding solutions to (see, e.g., [2]) 



(2.7) 



AJ = s{s--]f 



where 



f{r,0) 



iW 



$(r). 



(2.8) 



Note that we will have to set re = A; ± i and / = A ± i in appropriate places to make 
use of (2.7) and (2.8) in the Dirac problem. In geodesic polar coordinates one finds 



R 



dr — - Kth ^ 
shr 2 



(2.9) 
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= "1"^' {^^ + ^^"^ ^^^^ ^} ' ^^-^^^ 

and 

^>-=m{9r + -T-dr + -772-^0" - T^T-de + t— -j— ^ • 2.11 
[ shr sh^r 1 + cn r 1 + cn r J 

We make a substitution to deal with (2.7) and the side condition (2.8) [2,6]. Intro- 
duce the change of variables, 

$(r) = t'/2(l - t)^'/^F{t), (2.12) 

with 

chr — 1 IzP ,^ 

t = = ^ (2.13) 

chr + 1 V ^ 

in (2.7) using (2.11). One finds that F{t) satisfies the hypergeometric equation 

t(l - t)F" + (c - (a + 6 + l)t)F' -abF = (2.14) 

with 

a= — -K, b= — + K + l, c = l + l. (2.15) 

Note that it is precisely to simplify the description of these parameters that s was 
introduced in (2.5) and (2.6) above. For c ^ integer, two linearly independent 
solutions to the hypergeometric equation are 

F{a,b;c;t) and t^~''F{a - c+ l,b - c+ 1;2 - c;t), 

where F{a,b]C]t) is the hypergeometric function of Gauss. Putting this together 
with (2.8) and (2.12) we see that two linearly independent solutions to (2.7) satis- 
fying (2.8) are given by 



(2.16) 



V2{T,e,l,K,s) :=e^'^PL/2,«(r-), 

where ^ 

Pl,n{r) •= ^ T{l-X ^^"^ + n,a-K-l,l-l,t). 

It will also be useful at this point to record the solution of (2.7) of the form (2.8) 
which tends to as \z\ R (or r ^ 00). Define 
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and write 

Then v{r, 9, 1, hi, s) is a solution to (2.7) which tends to zero as r 
connection formula for the hypergeometric equation one finds 



(2.17) 

oo. From the 



v{l, K, s) 



4sin(7r/) 



V{a- n)T{a + K-l) 

V2{1,K,S) - -vi{l,K,s) 

V[a + K)V[a — K — I) 



where a = ^ and we have written v{l, k, s) for v{r, 6, /, s). 

Because we are interested in using the functions Vj and v to construct solutions 
to the Dirac equation, it will be useful to compute the action of the operators K,^ 
and K* on the functions Vj and v. Writing Vj{l, s) for the function Vj{r, 9, 1, k, s) 
we now summarize the action of Kk, and Kt on these functions: 



R 



R 



Vi{l - 1, K + 1, s), 

{^ + k){^-k-1)v2{1-1,k + 1, s), 



K,v{l, k,s) = ^{^ + k){^-k-1) v{l - 1, ^ + 1, s), 
K*vi{l, k + 1,s) = -^{^ + k) -K-l)vi{l + l,K, s), 

K*V2{1, K+1,S) = -^V2{1 + 1, S), 
K*v{l, K + l,s) = —^v{l + 1, K, s) . 



Finally we introduce the basic wave functions that we will use to analyze the solu- 
tions to the Dirac equation: 



vi{l + l,k- |,s) 

V2 {-I - |, /C + ^, s) 
^V2 {-l+^,k-l,s) 



(2.18) 



and 



Wlk 



( 



(-Z + |, /c - |, s) 



where s is the solution, 



s = R-^ (l + ^Jrn^R^ + 4fc2) , 

to (2.5). The solution wi^^ which is well behaved at infinity is 

r(a - /c + |)r(a + /c - 2 



(2.19) 



—ml 



Wl,k = 



4cos(7rZ) 



l,k 



T{a + k- ^)T{a -k + l)mR 
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W-l,k 



In order to see that these are indeed solutions to the Dirac equation it is useful 
to recompute the action of K^^ and K* on vj and v in the case where s is given by 
(2.5). Using 

2. (Rs _ 1 _|_ (Rs _ I — m'^R 

R \ 2 2^ '^J \ 2 2 '^J ~ 2 



one finds 



K,,_iv{-1 + i, A; - I) = ^v{-l - i A; + I), 



2 
2 

R 
2 



- i /c + i) = + I, /c - i) . 

We have dropped the explicit dependence on s, it being understood that s is given 
by (2.19). 

As mentioned in §1 the operators Mi and M2 are ladder operators on the 
eigenvectors for M3. Next we record the action of these infinitesimal symmetries of 
the Dirac operator on the wave functions wi^k ^ind wfj^: 

M^wi^k = {k + l)wi^k, 
Mswlk = {k-l)wlk, 

Miwi^k = wi-i,k, 
M2wlk = wl_i^k, 

M2Wi,k = m2{l, k)wi+i^k, 



(2.20) 



where 



mi{l,k) := (h^ + P-(/c-Z-|)'), 
m^il, k) := + e-(k + l+l)^y 



These results provide the rationale for the normalization of the we have chosen for 
the wave functions. The wave functions wi^k and are parametrized so that each 
becomes less singular at 2; = as Z increases. This is the reason for changing the 
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sign of / in the definition of Wi j^. The action of those infinitesimal symmetries that 
increase the local singularity of the wave functions will play an important role in 
what follows and the normalization we have adopted makes the action of Mi on wi^k 
and of M2 on w*^ as simple as possible. Finally we have chosen the normalization 
for wli^ so that the following conjugation symmetry is valid, 



where C 



This conjugation will play a role later on and it is especially 

important in the case k — where it is the conjugation with respect to which the 
Dirac operator is a real operator. 

Wc will now make use of the covariance of the Dirac operator under the action 
of SU(1, 1) to center these wave functions at the point a e D/j. Before we do 
this we will establish a convention for dealing with the branches of the multivalued 
functions wi^k and wf i^. Suppose that one has chosen a geodesic ray, £, at the point 
a = in Dj?. Then we fix a geodesic polar coordinate system with {9 = 0} — i 
and on Di^\£ we choose the branches of wi^k and w^*^ that are associated with the 
determinations (2.16) and (2.17) for vi and V2 with < ^ < 27r. 



Now consider the translation to the point a by the rigid motion T[a] on Dji given 
by 

T\a]z = = + . (2.21) 



-- + 1 

RR^ ^ 



az 



Then one easily computes that 



v{T[aiz) 



— az 



R^ 



az 



As a function of z this clearly does not wind around the origin and we can thus 
uniquely specify a branch of the logarithm by 

logw(T[a],0) = 

for all a G D^^. Using this choice of the logarithm to determine the fractional powers 
of V in the following formulas, we now define the local wave functions centered at a 

by 



v{T[-a],z)''- 



wi^k{z,a) := 



v{T[-a], z)''+^ 




wi,k{T[-a]z) 



(2.22) 



and 



v{T[-a],z) 






v{T[-a],z) 



<fe(T^[-«]^)- 



(2.23) 
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Note that a] occurs in this formula not T[a]. Because the functions wi^k and w^*^ 
are multivalued the formulas (2.22) and (2.23) do not make sense as they stand. 
We make sense of them in the following manner. Suppose that £ is a geodesic ray 
at a. Then T[—a]~'^£ = T[a]i is a geodesic ray at and hence determines a branch 
for wi^k and wfj^. With these branches for wi^k and wfj^ understood in (2.22) and 
(2.23) the functions wi^k{z, a) and i^{z, a) are well defined for z ^ £ and they are 
branched along £. We will think of wi^k{z, a) and Wi i^{z, a) as multivalued functions 
with specific branches determined by the additional specification of a geodesic ray 
at 2; = a. 

The importance of these wave functions is that the multivalued solutions to the 
Dirac equation branched at "a" that are of interest to us will have local expansions, 

W{z)= ^ {anWn+\{z,a) + bnW^_)^{z,a)} . 



Differentiating local expansions with respect to the branch points. In this 
subsection the dependence of the wave functions wi^k and wf on the parameter k 
will not be very important and so we will abbreviate wi^k = wi and wfj^ = wf. In 
order to diff'erentiate local expansions in the a parameter it will be useful to develop 

reciprocity relations for the behavior of wi{z, a) and wf{z, a) as functions of z and 
a. Let a(e) = T[eR]a and z{—e) = T[—eR]z. One easily computes that 

v{T[-a{e)],z) = ^Ii^^±^v{T[-a], z{-e)) 
rt — €Z rt -\- €(1 



and 



R + ea 

Observe also that under rotations wi and wf transform as follows, 



gi(fc+|)e Q 
e*^''-^)^ 



wi{e''z)^e 



i{k+i)e 



wi{z) 



and 



gi(fe+i)e 






gi(fe-i)e 



wt{e''z)^e 



■.{k-i)t 



wtiz). 



Of course, neither of these equations is correct without the proper interpretation of 
the multivalued functions wi and Wi . It will suffice for our purposes to note that the 
equations are correct provided that wi and Wi are branched along the geodesic ray 
£ at 0, we suppose that z ^ £, and we suppose that 9 is restricted to a sufficiently 
small neighborhood of e R so that the curve 9 ^ e^^ z does not intersect £. 
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Now we write 
wi{z,a{e)) = 



v{T[-a{e)], z)''+^ 

v{T[-a{e)],z)^-^ 



wi{T[-a{e)]z) 



and make use of the expressions for v{T[—a{e)], z) and T[—a{e)]z above. In the 
formula for wi{z, a(e)) that results, the transformation of wi under rotations implies 



Xk+h 
A'^-i 



wi {XT[-a]z{-e)) = X^+'wi{T[-a]z{-e)) 



where 



A 



R + ea 
R + ea 



Furthermore one has 

' v{T[-a], z{-e))''+'^ 




viT[-a],zi-e))'^-i 



wi{T[-a\z{-e)) = wi{z{-e),a). 



The reader should understand these formulas in the following manner. A geodesic 
ray £ at a has been chosen so that z ^ £ and e is chosen in a sufficiently small 
neighborhood, ?7, of e C so that the range of 

R + ea , . 

U 3e^ z{-e) 

R + ea ^ ' 

does not intersect £. The functions wii^z^a) and w\{z^o) are understood to be 
branched along £ and the functions wi and w\ are understood to be branched along 

T{-a\l. 

Putting the preceeding results together one gets the reciprocity formula. 



wi{z,a{e)) 



R + ea 



k->rl r 



R + ea 

where we made use of the fact that 

v{T[-eR],z) = 
In a precisely similar fashion one finds that 



v{T[-eR],zf+^ 

v{T[-eR],zf-h 



wi{T[-eR]z,a) 
(2.24) 



R — ez 
R — ez 



w!iz,a{e)) 



R + ea 
R + ea 



k-i 



v{T[-eR], z)''+^ 




v{T[-eR],z)''-^ 



wl{T[-eR]z,a). 

(2.25) 
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These formulae will allow us to compute the local expansions for the derivatives of 
the wave functions wi{z, a) and Wi{z, a) in the a variables. 

One reason for the utility of these reciprocity relations for differentiating the 
wave functions wi{z,a) and Wi{z, a) are the following easily confirmed formulae, 



and 



d_ 
d~e 



d_ 



v{T[eR],z)''+^ 

v{T[eR],z)''-^ 



v{T[eR],z 




v{T[eR], z)''-^ 



F{T[eR\z)U=o = MiF{z) (2.26) 



F{T[eR]z)U=o = M2F{z). (2.27) 



We will now use the reciprocity formulae above to provide local expansions for 
the derivatives of the wave functions wi{z, a) and Wi{z, a) in the "a" variable. The 
chain rule gives 



d^wi{z, a(e))|g=o = R ( daWi{z, a) - -^daWi(z, a) ) , (2.28) 



and 



deWi(z,a(e))\e=o = R ( daWi(z,a) - -^daWi(z,a) ) , (2.29) 



with completely analogous formulas for 'wf{z, a). Differentiating the right hand side 
of the reciprocity formulas (2.24) and (2.25) with the use of (2.26) and (2.27) one 
finds, 

d,wi{z,a{e))\,=o = {k + l)^wi{z,a)-MiWi{z,a), (2.30) 

Jri 

and 

— a 

d^wi{z, a(e))|e=o = -{k + l)—wi{z, a) - M2Wi{z, a), (2.31) 

Jri 



with analogous results for tyf (-2, a). To compute the action of Mi and M2 on wi{z, a) 

it is useful to represent Mj for j = 1, 2,3 in terms of the centered operators M^""^ 
obtained by conjugating Mj with the action (2.23) of r[— a]. One finds that, 

M2 = ^ ^ [uf - }?Mf^ - Ihuf^ } (2.32) 



M3 



= _J_ I^m}") - 6Mi"^ + (1 + |6|')Mi"^} , 
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where we have written 

for the scaled variable. Now substitute these results in (2.30) and (2.31) and equate 
the results with (2.28) and (2.29); one can solve for di,wi{z, a) and di,wi{z, a) to find 



and 



dbWi{z,a) = -^i^^^wi{z,a) - ^rq^ (^M^"''>wi{z,a) - ^ ^^^^^^ M^^^wiiz, a)^ 



The analogous results for Wi{z, a) are 



dbwiiz, a) = |^p ^r(^, a) - YZT^ yM^^'wliz, a) + j^p ^s ^H^^ «) 
and 

d,wt{z,a) = -^^^wtiz,a) - (^Mi^^wt{z,a) - Y^Mt^^t{z,a)^ . 

It is now a simple matter to use the results (2.21) for the wave functions to get the 
local expansions, 

dbwi{z,a) = - JZ^^'-i^^' 



{k + l)b m2(/,/c) . . 

dbwi{z,a) = -^—j^wi{z,a) - y— y^«;/+i(2;, a), 



o */ N {k-l)h mi(/, /c) ^ 

a6«;; (2;, a) = - ^^p W; (2;, a) - ^ _ W;+i(2:, a), 



(2.33) 



Estimates at infinity. We conclude this section with results that we will use to 
justify some later applications of Stokes' theorem in the non-compact domain T)r. 
If one considers (2.9) and (2.10) or (2.11) then one finds 

4 I ch. 1 xtv /sj^ 
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Next we fix s to the value given by (2.19). Then we introduce 

for brevity and also because $/,fe,s(^) is real valued. Then 

^;(r,^,Z,A;,s) = e^'(^+-)$«,,,,(r). (2.34) 

The differential equation (^K*K^ + m^) v = becomes the ordinary differential 
equation 



^iK s + -;— * Ks - { 9-^^ 9-^^ + + 1) H : — )^iKs = ^ 

1,1^,8 ^ l,K,s I ch^(0 Ch2(0 ^ ^ " ' ' 



for the radial wave function $i,K,s- We eliminate the first derivative term in this last 
equation with the substitution 

ui,K,s{r) = Vshr $/,«,s(r). 
After a little computation one finds that 

(2.35) 

As r — cxD the coefficients of m^^k.s in the third and fourth terms on the left hand 
side of (2.35) tend to exponentially fast. It is not surprising then that the behavior 
of the solutions to (2.35) is governed by the first two terms. When 

— \-k{k+1) + - > 

4 4 

there is one solution to (2.35) which is exponentially small at oo (in Hartmann's 
terminology [3] this is the principal solution). Suppose now that k = k — ^. Then 

= ^^ + k{k+1) + - = ^^ + k^ (2.36) 

is clearly positive and it follows from Theorem 8.1 in Hartmann [3] that the principal 
solution g to (2.35) has the asymptotics 

Wi,fe_i,s(r) = 0(e-^^) asr^oo. 
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Thus 

^^fe-i,.(0 = O(^^J asr^oo. (2.37) 

We are mostly interested in the consequences of (2.37) for the solutions to the 
Dirac equation which arc well behaved at oo. Since the second component of wi^k, 
for which we write {wi,k)2J J^^^ v{—l + ^, /c — |, s) it follows that 

/ -Mr\ 

I Mr))^ I = |$_^+i,fe_i,,(r)| = O as r ^ oo. (2.38) 

We would like to have something similar for the first component (wi^k).^. This 
component is a solution to the Helmholtz equation 

[K,_.K;_,+m')w = 0. (2.39) 

Since, 

S(k + -) 

it follows from (2.35) that the leading behavior (at oo) of the second order ODE 
that governs the radial part of the wave function Vshr (i«i,fc)i is 

u" - M'^u = 0. 

Then one has the following asymptotic estimate for the first component of wi^k-, 

mR ( e~^'^ \ 

|(w;/,fe(r))i| = ^|$_i_i,fe+i,,(r)| = I -^^== 1 asr^oo. (2.41) 

We will now use the estimates (2.38) and (2.41) to establish the following represen- 
tation for the norm, 

/ |tu/,fe(^)|^cZ//(^) = Q,fe(ro), (2.42) 

where 

Q,fe(r) = 27ri?2sh(r) i i ,,(r) 1 i ,,(r) (2.43) 

and we write 

Da,h = \z : i?tanh ^ < < -Rtanh - 
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To prove (2.42) suppose first that / and g are multivalued solutions to the Dirac 
equation in Dij\{0} with the same unitary monodromy about 0. Then a simple 
calculation shows that 

m{7 ■ g) dfJ^iz) = 2(i I -J^^idz 

Since f^gj is a smooth single valued function away from the origin Stokes' theorem 
implies 



m 
2" 



/ f-g {z)dii{z) = f -^^idz - I -^^idz, (2.44) 



where we write 

Cr = ^z: z = i?e*^tanh^|. 

In (2.44) let / = (7 = wi^k- Then substitute geodesic polar coordinates (2.0) in 
(2.44). Finally use the asymptotics (2.38) and (2.41) in the resulting equation to 
evaluate the limit ri 00. One obtains (2.42). The reader should note that the 
sign change due to the factor e*'^' by which k differs from « is important 

in getting the signs to work out properly in (2.42). 

Now suppose that in the exterior of the the disk 1^1 < — e the function / is 
a multivalued solution to the Dirac equation with monodromy multiplier e'^'^^^. If 
/ is square integrable in the exterior of |^| < tq — e with respect to the measure dfi 
then it has an expansion 

(2.45) 

n€Z+i 

This expansion gives the Fourier series representation for / in the angular variable 
9. If one makes use of this to do the circuit integrals that appear on the right hand 
side (2.44) then one finds, 



m 
2" 



[ J^idz= |cn(/)|'c„+A,fc(r). (2.46) 
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From (2.42) it is clear that Cn+x^k{'>^) is a positive, monotone decreasing function of 
r which tends to as r — > 00. Thus the monotone convergence theorem (for sums) 
and (2.46) imply that 

lim f -ll^idz = (2.47) 
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when / is a multivalued solution to the Dirac equation which is square integrable 
in a neighborhood of infinity. Finally, we want to know that something similar is 
true if / and g are two multivalued solutions to the Dirac equation both square 
integrable in a neighborhood of oo. Consider the Hermitian form, 



J?2 

Then (2.45) shows that this is positive definite and so the Schwarz inequality, 

|(/,^)rP < U,f)r{g,g)r, 

is valid. It follows immediately from (2.47) that 

lim I -^^^idz^O. (2.48) 

We summarize these developments in the following proposition: 

Proposition 2.0 Suppose that f and g are two solutions to the Dirac equation 
which have expansions of type (2.45) in the region Dro.oo- Then 

lim / •^^^^.o idz = 0. 



Jr 1- ^ 



§3 Existence Results 

A model for the simply connected covering Di?(a). To understand the asymp- 

totics of the Green function for a Dirac operator with vertex insertions it will be 
important to have an understanding of certain spaces of multivalued solutions to 
the Dirac equation. We will follow SMJ [15] here in establishing an existence 
theory for wave functions with specified branching and restricted singularities. Let 
{ai, 02, ... , On} denote n distinct points in the hyperbolic disk D^. It will be useful 
for us to use complex variables writing aj = aj + iPj when aj = [aj, We will 
also write 

a = {tti, a2, . . . , cin}- 

Following SMJ [15] we will now show that the space of solutions, w, to the Dirac 
equation, 

(m - Dk)w = 0, 



28 



which are branched at {ai, 02, . . . , On} with monodromy e^^*'^-' at aj (with < 
|Aj| < ^ ) and globally in L^(Dij, djj) is an n dimensional space. 

It will be convenient (particularly in discussing smoothness questions) to work 
in the simply connected covering space of Dfj\{ai, 02, . . . , fin} and we will now spec- 
ify our conventions for such considerations. Write Dfl(a) for Dfl;\{ai, a2, . . . , a^} 
and Di?(a) for the simply connected covering space of the punctured disk Di?(a). 
Multivalued functions on Djj(a) can, of course, be regarded as functions on Djj(a). 
Introducing appropriate branch cuts, at aj a multivalued function is also spec- 
ified by a particular branch defined on Djj\{£i,£2) • • • )^n}- Both descriptions will 
be important for us and to pass easily from one to the other we will always work in a 
specific model for Y)^{a). Fix a base point ao in Djj different from any of the points 
ai, a2, . . . , an- We take our model of Di^(o) to be the homotopy classes of paths in 
Dij(a) which start at oq. The projection, pr, from Dij(a) onto Djj(a) maps a path 
7 onto its endpoint 7(1). This model of Dit(a) comes with a distinguished point, 
ao, in the fiber over oq which is the class of the constant path t — > oq. 

For simplicity we fix the convention that our branch cuts Ij are always geodesic 
rays joining the point aj to a point on the boundary of D^^. If {£1, £2, • • • 5 ^n} is a 
pairwise disjoint collection of branch cuts then the inverse image of 

Di?\{£i,£2,...,M 

under the projection, 

pr : Di?(o) Dfl, 

splits into path components on which the projection is a diffeomorphism. Let Co(£) 
denote the path compojicnt of this inverse image which contains the point Sq. If 
/ is a map defined on Dij(a) then we will refer to the restriction of / to the path 
component Co(^) as the principal branch of /. The principal branch of such a 
function / can also be regarded as a function defined on Y)f{\{li,l2i ■ ■ ■ J-n} and 
we will do this without further comment. Conversely if / is a function defined 
on Dfl\{£i, £2, • • • 7 -^n} we can regard it as a function on Co(-^). For multivalued 
functions that have a natural extension to maps on Dij(a) this gives a unique way 
to go from a branch for such a function to a function defined on Dji{a) . For example, 
we have already described how to fix a branch of the multivalued function wi{z, aj) 
when we are given a branch cut, £j, at aj. We can regard this as a function on Co{£) 
which has an obvious extension to a smooth function defined on Dij(a). When we 
wish to regard wi{-,aj) as a function on T)R{a) we will write wi{z,aj) with the 
understanding that is a point in Dij(a) which projects onto z e Di?(a). It is 
understood that a branch cut £j must be chosen before the function wi{z,aj) is 
uniquely determined. 
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Let [aj\ denote the homotopy class of the simple ao based loop, cij, about aj 



in Dij(a) which circles in a counter-clockwise fashion and does not wind around 
any of the other points with i ^ j. Let Rj denote the deck transformation on 
Dfl;(a) which maps ao to the point in Dfi{a) which corresponds to the homotopy 
class [ctj]- To simplify notation in what follows we will, as above, write z for points 
in Dfl;(a) and z = pr{z) for their projections in Dij(a). Suppose that / is a smooth 
function on Djj(a) for which 



If F is the principal branch of such a function / for some choice of branch cuts 
{£±,£2, ■ ■ ■ , £n} with ao ^ £j for J = 1, 2, . . . , n then we will say that F is branched 



at Uj with monodromy multiplier e 



2TviXi 



for j = 1,2, . . .,n. 



Multivalued solutions with specified branching. We now turn to the consider- 
ation of multivalued solutions, w, to the Dirac equation with specified monodromy 
at the branch points aj. We can locally expand such functions near the branch 
points in Fourier series in the angular variable 9 of geodesic polar coordinates suit- 
ably modified to reflect the appropriate monodromy. As in SMJ [15] these Fourier 
series translate into local expansions 

w{z) — al^{w)wri+\^{z^aj) +hi^{w)w'^_^^.{z^aj) lor z neox aj., (3-0) 

nez+i 

where (3.0) is given a precise meaning when the choice of a pairwise disjoint collec- 
tion of branch cuts {£±,£2, ■ ■ ■ ,£n} has been made with ao ^ for j = 1, 2, . . . , n. 

The following local asymptotics for the wave functions may be deduced from 
the formulas for wi (z) and Wi {z) in the preceeding section: 



wi{z) = 



2e^ 



4). 



(§)' 



mi?r(Z+i) 



+ 0(r 



(5) 



1+4 



r(«+T) 



+ 0{r^+i) 



(3.1) 



and 



2e" 



+ 0(r'+2) 
- + 0(r'+^) 



(3.2) 



mRr{l+^) 

From these asymptotics we deduce that wi and are locally in L'^idfj) provided 
that Z > — |. Thus in the expansion (3.0) of an wave function we expect that 
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the coefficients a^(w) will vanish unless n + Aj > — | and the coefficient Vn{w) will 
vanish unless n — Xj > —\. This leads us to define 

f a\ (w) if Xj < 
cAw) ^ { . ^ . ^ (3.3) 



and 



if > 



The coefficients Cj{w) and c*(t(;) are the lowest order coefficients that can occur 
(i.e., be non zero) in an expansion of type (3.0) for a function which is locally in 
{djj) . We will now make this more precise with a formula for the inner product 
of two multivalued solutions to the Dirac equation with the same monodromy. The 
infinitesimal version of formal skew symmetry for is 

fl92 .,_ f29l 



(/ • Dkg + Dkf ■ g) d^i{z) = 2d < ^^^^ \^^^^ \ ■ (3-^) 

The d on the right hand side is the usual exterior derivative. This is true if / and 
g are smooth functions on D^j and it will remain true for functions / and g that 
are smooth branched functions with the same unitary monodromy at the points aj 
for J = 1, 2, . . . , n — at least if one stays away from the branch points. In this case 
each side of (3.5) will be single valued on Di?\{ai, 02, . . . , a^}. 

Theorem 3.0 Suppose that f and g are multivalued solutions to the Dirac equation, 

(m - Dk)f = (m - Dk)g = 0, 

which are branched at ai, 02, . . . , with monodromy 6^"^*^^ at aj. If f and g are 
in L?'{T>ii,dji) then 



rrr'R 



/n n 
-^R 7 = 1 7 = 1 



where Sj = sinTrAj. 

Proof. Substituting D^f — mf and D^g — mg in equation (3.5) one finds 

'm(f-g) dii{z) = d \ -^^^^id:z ^^^idz 
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In fact, calculating the exterior derivatives of each of the terms on the right hand 
side of this last equation one finds that more is true in case / and g satisfy the 
Dirac equation, namely, 

^{7- g) dii{z) = 2(i I ^ ^^%\2 ^^^1 = -2<^ | ^ '^^^^p ^dz | . (3.7) 

Now we wish to calculate the inner product of / and g by first integrating 
(3.7) over the complement of the union of the disks of radius e about aj for j — 
1, . . . , n inside a disk of radius p. Then we use Stokes' theorem to reduce the integrals 
on the right to the sum of integrals over the circles, Cj(e), of radius e about the 
points ttj with the standard counter-clockwise orientation (this will account for 
some minus signs showing up below) and an integral over the circle of radius p. It 
follows from Proposition 2.0 that the boundary integral over the circle Cp is in the 
limit p — > oo. Finally we recover the inner product from the remaining circuit 
integrals in the limit e ^ 0. The result is 

^ I f-gdp{z)=— lim j -^^^^idz 

Now introduce the local variables Uj — T[—aj\z in the Cj{e) integrals. The change 
of variables gives 



1 - ^ 

idz = —, . idu. 



and 



1 1 ""j 



idz — — ^ idu 



'3 



fl2 



•J- 



Substituting these results in the appropriate contour integrals above and taking 
limits as e ^ of some of the multiplicative factors that appear inside the integrands 
one finds, 



m 
~2 



/ f-9dpiz)=- lim y2 fi92iuj)iduj 



3 

n „ 

= lim V / f2gi{uj)idu 
7=1 



J' 

e ^ j^l-^Cjie) 
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where the notation fjQkiu) is meant to remind the reader that it is the product 
fjgk that is a function on Dij(a). Geodesic polar coordinates {rj,9j) about aj 
(determined by the branch cut £j) give 

^ 2 2 

We can use this asymptotic formula in the contour integrals to get 



m 
2" 



/ f-g dfi{z) = - hm hg2{uj)e-^'^dej 
e ^ 7=1 ^ 



Finally if we use the asymptotics of the local wave functions, given by (3.1) and 
(3.2), in these last two formulas and the identity, 

^ 7 

^ = cos TTt , 

r(i-or(i + 

then one obtains the two local expansion results for the I? inner product given 
in (3.6) above. To see this it is also useful to note that the multiplicative factor 
f (T[— aj], 2;), which appears in the wave functions, wi{z,aj) and wf{z,aj) tends to 
1 as z ^ aj. QED 

Definition. Let I denote the n x n diagonal matrix with jj entry Ij with Ij chosen 
so that < Ij < 1. Define W{a,l) to he the space of multivalued solutions, w, 
to the Dirac equation (m — Dk)w — which are branched at {ai, 02, ... , an} with 
monodromy 6^"^*'^ in a counter-clockwise circuit of aj and which are globally in 

Remark: We parametrize the monodromy in this definition by Ij chosen between 
and 1 rather than Xj chosen between — | and | since the local expansions of func- 
tions in W(a, I) will be slightly simpler to characterize in terms of the parameters 
Ij. 

Theorem 3.0 shows that the elements, w, in W(a, /) are completely determined 
by the coefficients Cj{w) for j = 1,2, ...,n or by the coefficients c*(io) for j — 
1,2, ...,n. Thus the dimension of VV(a, /) is less than or equal to n. Following 
SMJ [15] we will now give a functional analytic proof that the dimension of >V(a, /) 
is exactly n by constructing a canonical basis. The functional analytic part of 
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the proof requires that we first consider an associated problem for the Helmholtz 
operator + K/^K* where n = k — ^. Recall from (2.2) that 



(m-Dfe)*(m-Dfe) = 



+ K^K* 
m2 + K*K, 



As a start towards constructing solutions to the Dirac equation with prescribed 
branching at {oi, 02, • . . , a^} we will first consider the problem of finding multival- 
ued solutions, V, to the Helmholtz equation 



(3i 



We will construct solutions of the desired type by fixing the leading singularity of 
V at one branch point. Then we will get a solution to the Dirac equation of the 
appropriate type by considering 



w 



ill ^K*_^' 



(3.9) 



The reader might observe in what follows that using (3.8) instead of 

(m^ + K:K,)v = 0, 

and (3.9) instead of 



w 



—m ^Kj,_iv 



is dictated by the desire to produce canonical solutions satifying Cj{Wi) = 6ij rather 
than the dual canonical basis c*{Wi) = 5ij. 

To formulate local expansions for the Helmholtz equation in a convenient form 
we wish to introduce here the analogue of the wave functions wi{z, aj) and Wi{z, aj) 
for Vj{l, K, s) introduced in (2.16) above. Recall from §2 that (3.8) is equivalent to, 



with K = k — h and s given by (2.5). Now we write 



and 



VlA^) = Vi{z,l,K,s), 
vlniz)=V2iz, -1,K,S), 



(3.10) 
(3.11) 
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with s given by (2.5). Then vi^^+i and Vi are both multivalued solutions to 

(3.8) above with monodromy e^'^*' and e"^'^*' respectively, in a counter-clockwise 
circuit of 0. We translate vi^^ and v*^ so that they are "centered at a e D^^" in the 
following manner, 

vi,^{z, a) := v{T[-a],z)''vi,^{T[-a]z) (3.12) 

and 

vtjz,a) := v{T[-a],zrvtjT[-a]z). (3.13) 

Precisely as in §2 a branch cut £ at a is used to fix a unique branch for the functions 
vi^k{z, a) and Vij^{zj a). The local singularity in both vi^i^{z) and ^^{z) goes like r' 
where r is the hyperbolic distance between and z. To make use of the formula 

(3.9) to get an solution to the Dirac equation we will require that the solutions 
we consider for (3.8) have both v and K*v in L^. It is now a simple matter to 
check (using the formulas for the action of K* on vi and v* given in the preceeding 
section) that if is a multivalued solution to (3.8) with monodromy multiplier e^'^*^^ 
in a counter-clockwise circuit of aj and if v and K*v are in then provided Ij is 
chosen so that < Ij < 1 the function v will have restricted local expansions 



v{z) 



CO 



oo 

ai{v)vn+i^,n+i{z,aj) + J2Pii'"K 

n=l 



(3.14) 



for z near aj. In each case the sum is now over integer values of n rather than half 
integers. Our strategy in finding such solutions will be to seek a canonical basis 
{Vi} with 

a^^{Vi) = 5ij. (3.15) 



Existence for a canonical basis. Let (fj{z) be a Cq° function which is 
identically 1 in a neighborhood of the point z = aj and which vanishes in the 
complement of an open ball about z = aj. This open ball should be small enough 
so that (fj vanishes in a neighborhood of each of the points with i ^ j. Observe 
next that if Vj is a solution to (3.8) with local expansions of type (3.14) and which 
further satisfies (3.15) then, 

(m^ + K^K*) {Vj - cpjvi.-i,^+i{z, aj)) = 
- (m^ + K^K*) (pjvi.-i,^+i{z,aj) :=fj, 

where fj is a smooth multivalued function which vanishes in a neighborhood of 
the point aj (and in a neighborhood of all the other branch points as well). For 
simplicity we will refer to branched functions that have smooth extensions (from 
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some Co{£)) to functions on Di?(a) as smooth multivalued functions. We will con- 
struct Vj by finding a suitably regular solution to (3.16) — one whose local expansion 
coefficients q;!_i vanish for z = 1, 2, . . . , n and which is in the Sobolev space H^. 
Motivated by (3.16) we are interested in the existence question for solutions to 

{m' + K^K:)uj^f, (3.17) 

for Uj in a completion of 

Co°°(Dfl(a)) 

with specified branching at for z = 1, 2, . . . , n. Here CQ°(Dii(a)) denotes the space 
of functions on Di^(a) whose supports project onto compact subsets of Dij(a). 
In particular such functions vanish in a neighborhood of the points ttj and in a 
neighborhood of the boundary of Dr. Let denote the subset of Co°(Djj(a)) 
consisting of those functions / which transform under the deck transformation Rj 

by, 

/(i?,i) = e2-'V(i), 
where < Ij < 1 for j = 1,2, . . . ,n and we write 

I — {h, h, ■ ■ ■ , In)- 

Now let Ha,i denote the Hilbert space completion of with respect to the norm 
derived from the inner product, 

{F, G)= [ {K*F(i)K:G{I) + m'F^G{I)] di,{z). 

Note that the integrand descends to a function on Dj^ because the monodromy 
multipliers of F and G cancel. Now consider the element /-,■ e Ca,i defined above 
in (3.16). For v e consider the linear functional, 

[ W)<^ M^) =< fj, V > . (3.18) 

Since {h, h) > w? < h,h >, for all h e C^^, it is clear that the linear functional 
(3.18) is continuous on the Hilbert space Ha,i- Thus by the Riesz representation 
theorem there exists an element Uj e Ha,i so that 

Since the formal adjoint of is K,^ one sees from this last equation that Uj is 
a distribution solution to (3.17). We wish to be a little more explicit about the 
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local regularity of this solution near the branch point points aj for j = 1, . . . , n and 
for this purpose it is useful to consider in more detail functions F e i^a,/- Such a 
function F, is the limit a sequence of functions e in the Ha^i norm. Recall 
that functions in have "compact supports" so that integration by parts in the 
formula for the Ha,i norm of Fm is justified and one finds that (F^, Fm) is, 

/ F\;^K,K:Fm{T, di^{z) [ FJ^F^{z) di^{z). (3.19) 

Now a simple computation shows that 

- = • 

Substituting — ^'^'^^^^ for Ki^K* in (3.19) and integrating by parts one 

finds that the L'^{dji) norm of K^^iFm can be expressed in terms of the L'^{dji) 
norm of K*Fm and the L'^{dfx) norm of Fm- The analogous result for the difi'erence 
Fm Fm' shows that K i^^iFm is a Cauchy sequence in L'^idji). Thus we find that 
Kft+iF is in L'^{diJ,). 

Next observe that if F, K*F, and K^+iF are in Lp'{dii) it follows that F is 
locally in the Sobolev space (at least away from the branch points). Returning 
now to the distribution solution Uj e Ha^i of (3.17) above one can now conclude 
from the standard regularity results for solutions to elliptic equations that since Uj 
is locally in the Sobolev space it must actually be C°° away from the branch 
points for z = 1, 2, . . . , n. 

Using the fact that tt^, K'^Uj and K^+ittj must all be locally in it is not 
difficult to see that the local expansions for Uj must have the form 

oo oo 
= XI «nK-)^n+«i(^, Oi) + (3l^{Uj)v*^_i.{z, a^) 
n=0 n=l 

for z near to Oj. The reader should note that it is precisely the additional condition 
that K^j^iUj must be in that forces a^_i{uj) = for all i. Thus if we define, 

= Uj{z) + >^j{z)vi.-i^^+i{z, aj), 

then it is clear that Vj is a solution to (3.9) which is appropriately branched at 
ai, a2, . . . , fin, and it is in L^(Dij, dfi) together with K*Vj by construction, since 
K*vi.-i is locally less singular than vi.-i. Finally 

a!_i(^,) = % 
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since Uj does not contribute to the local expansion coefficients C(-i{Vj). To exhibit 
the dependence of Vj on the various parameters we will write, 



where I = {li,l2, ■ ■ ■ Jn)- 

We now exhibit the canonical basis of Dirac wave functions. Define, 



Wj{l, k) 



2 ^ ^ 



(3.20) 



where < Ij < 1 for each j = 1, 2, . . . , n. Then (i, k) is an wave function for 
the Dirac equation with monodromy e^'^''^ about the point aj and with 



Ci{wj{l, k)) 



'13 ■ 



(3.21) 



The response functions. The canonical wave functions are very useful but they 
are not the wave functions that are fundamental in our treatment of r— functions. 
We turn now to the consideration of the Dirac wave functions Wj(J, A, fc) and 
W* {z,\,k) that will be central for us. Because of the role they play in the for- 
mula (4.50) for the derivative of the Green function we will refer to these wave 
functions as the response functions. We start with a characterization of the desired 
response function Wj and then we use the wave functions to prove the existence 
oiWj and W* . 

The response function Wj whose existence we wish to establish is characterized 
by the following three conditions: 

(1) Wj{z, A, k) is a multivalued solution to the Dirac equation which is in L^{dfi) 
in the complement of some compact neighborhood of {ai, . . . , a^}. 

(2) Wj{z, A, k) has monodromy multiplier e^'^*^^ in a counter-clockwise circuit of aj 
(with Xj now chosen so that < \Xj\ < ^). 

(3) Wj {z, A, k) has local expansions at each of the points aj of the following type. 



Wj{z, A, k) = i+A, + Yl + Kji\ k)wl_^^ ] , (3.22) 

n>0 

where n G Z + i and we abbreviate Wn+\^^k{z, Ui) by Wn+Xi and w^_^, j^{z, ai) by 
W^_^.. The first thing to observe about (3.22) is that it characterizes Wj{z, A, k). 
Any two functions with local expansions of type (3.22) differ by a wave function 
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whose local expansion has only terms involving Wn+Ai and 'W^_^. for n > |. A little 
computation using Theorem 3.0 shows that any such wave function, w, has norm 
since either Ci{w) = or c*{w) = for each i = 1, 2, . . . , n. 

Remark. Because Wj (z, A, k) depends on many parameters we will often omit 
those which are of less interest in the current calculations and rely on the charac- 
teristic names for the variables to save the reader from confusion. Thus in what 
follows we will write Wj{X, k) for the function Wj{z, A, k). 

To prove the existence of Wj{X, k) we will consider two cases. Let P denote the 
subset of {1, 2, . . . , n} which consists of those j for which < Aj < ^ and let P' 
denote the complementary subset of {1,2, ...,n} consisting of those j for which 
— i < Aj < 0. If j e P then it is not hard to see that any function Wj{\, k) with 
local expansions of type (3.22) is locally in L^. In fact, one can see that if Wj{\ k) 
satisfies (3.22) and j G P then 

c^{Wj{\ k)) = 5ij for e P (3.23) 

and 

c*{Wj{X, /c)) = for i e P'. (3.24) 

But Theorem 3.0 shows that if -u; G W(a, I) and Ci{w) = for z G P and c*{w) = 
for z G P' then w = 0. Since the dimension of W(a, /) is n it follows that {ciji^p U 
{c*}j£P/ is a collection of independent linear functionals on W(a, I). It follows that 
there exists an element Wj(A, k) of Vy(a, I) which satisfies the conditions (3.23) and 
(3.24). Without difficulty one can see that this implies Wj (A, k) has local expansions 
of type (3.22) and finishes the existence proof when j G P. 

When j G P' the leading singularity in (3.22) at aj is no longer locally in L^. 
To prove the existence of Wj{\, k) in this case it is probably simplest to return to 
the construction of the function Vj . Instead of forcing the singularity of Vj to agree 
with vij-i^n+i{z,aj) near z = aj one can construct an analogous function Vj whose 
singularity at 2; = aj is matched with vxj-i^^+i{z, aj) where Xj = Ij — 1. From Vj 

we construct _ 

" ^V- 1 

as before. Again using the independence of the collection of Ci for z G P and the 
collection of c* for z G P' we see that there exists an element Gj G >V(a, /) so that 

Ci(Gj) = Ci{Wj) for z G P 

and ^ 

c*{Gj) = c*{Wj)iorieV'. 



Wj = 
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If we now define 



Wj{X, k) = Wj - Gj for K = k 



1 

2 ' 



then it is easy to check that Wj{X,k) has local expansions of type (3.22). This 
finishes our existence proof for the response functions satisfying (3.22). 

The response functions Wj{z,\,k) are characterized by conditions precisely 
analogous to 1-3 above. In particular W* is a multivalued solution to the Dirac 
equation with monodromy multiplier e'^'^^^^ in a counter-clockwise circuit about aj. 
Instead of (3.22), however, the local expansion for W* is. 



W;{z, A, k) = Sijw*_^_^^ + J2 + ^X-aJ • (3.25) 



The proof of the existence of the response functions W* is precisely analogous to 
the existence proof for Wj and we leave this to the reader. The definition of Wj we 
have given here differs by the change A — > —A compared to the analogous response 
function in [12]. In the present context the response functions Wj and Wj are both 
solutions to the same Dirac equation with the same monodromy. 

Remark. We would like to point out one possible source of confusion for readers of 
this paper and the original papers of SMJ. In [12] there are two different bases for 
the same space W{a, I) referred to as the canonical basis wj and the dual canonical 
basis w* . The response functions W* above should not be confused with a dual basis 
of any sort. As the reader can easily see, the span of the response functions Wj and 
Wj for j — 1,2, ... ,n is 2n dimensional as the local expansions (3.22) and (3.25) 
show that these functions are linearly independent. 

§4 The Green function G"'^ for the Hyperbolic Dirac operator 

The goal of this section is to construct a Green function for the Dirac operator in 
the presence of branch type singularities at the points aj for j = 1, 2, . . . , n. Before 
we attempt this there are a number of preliminary considerations. 

Green functions in the absence of branch points. We begin by identifying 
the Green functions for the Helmholtz and Dirac operators in the absence of any 
singularities. To start we require Green's identity. First observe that the formula 
which identifies K* with the formal adjoint of K,^ may be written, 



n>0 



{KJ{z)g{z) - f{z)Klg{z))d^i{z) = d<p{f,g), 



(4.1) 



where the one-form (p is, 




(4.2) 
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and d denotes the exterior derivative. Taking complex conjugates of (4.1) one finds 

(K*fg - fK^g)dfi = -d^^- (4-3) 
From (4.1) and (4.3) one may deduce 

{^K,g - jK:K,g)dfi = d<f(f, K^g) (4.4) 

and 

{KWJa - KjK^g)d,i = -dip{KJ,g). (4.5) 
Adding (4.4) and (4.5) we find the infinitesimal version of Green's theorem 

(AJg - 'fA,g)dfi = di<f{KJ,g) - <fi(f, K^g)) , (4.6) 

where we recall (2.3). Finally if we write, 

/ 2\ 2 4^«(k + 1) 

then we can rewrite (4.6) as 

{{A,-X{s))fg - 7(A« - X{s))g)df, = d{ip{KJ,g) - ^{J, K,g)). (4.7) 

We are interested in the Green function for A(s) — A^^ and so we will look for radially 
symmetric solutions to 

(A(s) - A,)^ = 

which have the appropriate singularity at 2; = and which are in Lp' (dfi) in a neigh- 
borhood of |2;| = R. Recall now from (2.7)-(2.15) that a solution to the Helmholtz 
equation can be written in the form, 

2 ~ 2 
9nA^) = (1 - ^) ' (^) ' (4-8) 

where -FK,s(t) satisfies the hypergeometric equation (2.14) with, 

and where s is again given by (2.19). 

The case c = 1 is special for the hypergeometric equation. One solution, 
F{a, b; 1; t) is regular at t = 0, while the second independent solution has a logarth- 
mic singularity at t = 0. A local analysis near t = 1 shows that when Rs = a + b > 1 
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only one solution, F^^s{t), determines a solution, g^^^siz), via (4.8) which is square 
integrable near \z\ = R (square integrable with respect to dfx that is). Of course 
this solution is only determined up to a constant multiple by this requirement. We 
fix the normalization of F^^sit) by requiring that 

g.A^) ~ log + 0{\z\ log 1^1) (4.9) 

for z near 0. For this to work we must know that the solution Fi^^s{t) which produces 
the solution g^^siz) which is square integrable near \z\ = R is not also the solution 
of the hypergeometric equation which is analytic at t = 0. We will sketch how one 
can see this. Suppose that v is the radial solution of 

which tends to as 1^1 ^ i?, then for k = /c — | 

'-K^_iv 



w 



2 

mv 



is a solution to the Dirac equation which is square integrable at oo. If v is also well 
behaved near \z\ = 0, then w will be square integrable in a disk, 

Dp = (z: \z\ < Rth^ 



2 



of radius p > about 0. A calculation similar to the calculation that produced 
(2.42) and (2.43) shows that the norm of w on Dp is expressible in a form 
precisely analogous to (2.43). The estimates (2.38) and (2.41) for w show that this 
expression for the norm tends to as p — > oo. This contradiction shows that the 
solutions well behaved at oo and at cannot match up. 

We next consider the sense in which g^^g is a Green function for A(s) — A^^. Let 
/ denote a smooth function of compact support in Dji. Integrate both sides of (4.7) 
with f = f and g = gK,s over the complement of a disk of radius e about z = 0. 
Use Stokes' theorem to reduce the result to a boundary integral over Ce, the circle 
of radius e about 0, and then use the asymptotics (4.9) to evaluate the boundary 
integrals in the limit e — > 0. We find 

/ iXis)-A^)f{z)'^Zi^dfi{z) = f{0). (4.11) 

This shows that gK,,s{z) is the solution appropriate to determine the Green func- 
tion for A(s) — A^. To obtain this Green function we must translate gK,s in the 
appropriate fashion. Introduce the variable 

_r n „ Rz — Ra 
w := T[-a]z = R^ 



-az 



+ R^' 
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and define 

, , R"^ + aw 
v{w) := „ _ ■ 
+ aw 

Let v{w)'^ denote the fractional power of v{w) normalized so that v{0)'^ — 1. Then 
the transformation property of which is the analogue of the transformation 
property (1-18) of the Dirac operator is 

^(^)-A^-)^H'^ = -K:{w)K^{w) + ^^^^^i^ := A^{w). (4.12) 

This should be understood in the following fashion. The operator A^™-* is the 
expression for A^ that one obtains by the change of variables z = z{w). The 
operators K^{w) and K*{w) are the analogues of and K* in the w coordinates, 



and 



R^ J i?2 

so that A^{w) is the same operator in the tu— coordinates that A^ is in the z 
coordinates. Now we introduce 

g,,,{y, z) = v{T[-y]zYg^,,{T[-y]z) = g.AA-v]^)- (4-13) 

To see that this is the appropriate Green function suppose that / is a smooth 
function of compact support in Djj as above and consider the integral 

/ 



(A(s) - A^)f(z)g^^s{y, z) dn{z) 

{X{s) - A^:'^)f{T[y]w)v{wYg,,s{w)d^i{w) 
D« (4.14) 



(A(s) - A,{w))v{w)-^f{T[y]w)g^,s{w) dii{w) 
= v{w)--f{T[y]wY=o = f{y), 



where we made the substitution z ^ w to obtain the transformed integral in the 
second line then used (4.13), (4.12) and finally (4.11) to obtain the third and fourth 
lines. 
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Next we construct the Green function for the Dirac operator, m — -D^, using 
the Green function gK,s{y-i z). Recall now (2.2) and (2.6) from which it follows that 



{m + Dk){m - Dk) 



X{s) 









A(^)-A,_i 



(4.15) 



Now let f{z) denote a smooth function of compact support on Tin with values in 
C^. Then from (4.15) and (4.14) we have 



/ {m + Dk){m-Dk)f{z) 







df^iz) = fi{y)- 



(4.16) 



If we exclude an e disk about z = y in (4.16) and use Stokes' theorem to "integrate 
by parts" once in (4.16) then we find 



/ 



(m - Dk)f{z) ■ Gi{y, z) dii{z) = /i(y), 



(4.17) 



where 



Gi{y,z) = {m + Dl) 







= (m - Dk) 







(4.18) 



The boundary of the e disk about z = y makes a contribution to the Stokes' theorem 
calculation that vanishes in the limit e — due to the weak logarithmic singularity 
in g^+i^siVj z) sX z = y. In a precisely similar fashion one finds that for 



G2iy,z) := (m - Dk) 







one has 



/ (m - Dk)f{z) ■G2{y,z)df,{z) = h{y). 



(4.19) 



(4.20) 



The Green function for the Helmholtz operator in the presence of branch 
points. Before tackling the existence question for a Green function associated with 
the branched version of m — D we will first consider the analogous problem for the 
"Helmholtz" operator 

X{s) - A« = + K^K,, 

with 

X / \ 2 , 4k(k + 1) 

A(g) = m^+ 
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We wish to consider the action of A„ on functions on the unit disk with specified 
branching at a finite collection of points {oi, 02, ... , a„} in the disk D^j. It will be 
useful for us to recall the conventions of §3. Thus we identify points aj = {aj,Pj) 
in the disk with their standard representation as complex numbers aj = aj + i^j. 
Write Dij(a) for Di?\{ai, 02, ... , a„} and Di?(a) for the simply connected covering 
space of the punctured disk T)ji{a). Let (7^(0^; (a)) denote the smooth functions 
on Dij(a) whose supports project onto compact subsets of Djj(a). Fix a base point 
ao in different from any of the points ai,a2, ■ ■ ■ ,an and let [aj] denote the 
homotopy class of the simple ao based loop, aj, about aj in Dit(a) which circles aj 
in a counterclockwise fashion and does not wind around any of the other points 
with i ^ j. Points in D/j(a) over ao correspond to homotopy classes of ao based 
loops in Dij(a) and as above we fix our model of Dij(a) to be the collection of 
homotopy classes of paths in Dji{a) which start at ao. The canonical projection, 

pr : Dfl(a) Br, 

maps each class [7] of paths that start at ao into the endpoint 7(1). This model of 
Dij(a) comes with a distinguished base point. Let ao denote the point in Di?(a) 
which corresponds to the constant path [0, 1] 3 t ^ gq. Let Rj denote the deck 
transformation on Dij(a) which maps So to the point in Dij(a) which corresponds to 
the homotopy class [aj]. Although it is not completely precise it will be convenient 
to adopt the following notation to avoid clumsy expressions. We will write z for a 
typical point in Dij(a) and 

pr(i) = z 

for the projection on D;^. 

We wish to consider the action of on a subset of C^(Dij(a)) with speci- 
fied branching at aj for j — 1,2, ... ,n. Let denote the subset of CQ^(Dij(a)) 
consisting of those functions / which transform under Rj by 

/(i?,i) = e-2-^'^7(i), 

where < /j < 1 for j = 1,2, ... ,n. Note that since is a module for the 
multiplicative action of smooth functions on Dr and since the differential operators 
dg and dz lift in an obvious way to differential operators on there is no difficulty 
in defining the action of on C^^. We will now construct a Green function for 
X{s) — acting on C^^. This is not completely accurate since we will implicitly 
limit the strength of the singularities at the points aj in a completion of but it 
will prove simpler to construct the Green function than it will to give a complete 
discussion of the domain of the corresponding operator. By analogy with (4.14) 
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above we seek a function g""' {y, i) so that for any / G we have 

/ (A(.)-A.)/(%«.Hy,i)dM^)=/(y). (4.21) 

Of course, for this integral to make sense it is necessary for the integrand, 

D^(a) 3z^ (A(s) - A«)/(%-.^(y,i), 

to descend to a function on JDr. This wiU happen if g"''''{y, i) transforms as foUows, 

g'''\y,Rjr, = e-^^''^g-\a,^, (4.22) 

under the deck transformations Rj. Equation (4.21) suggests that g°''''{y,z) should 
transform in the first variable 

g^'\Rjy,T, = e^-''^g-'\y,r,. (4.23) 

We will now give a functional analytic proof for the existence of a suitable 
version of g"''''{y, i). To begin we will first concentrate on the behavior of g"''\y, z) 
as a function of the second variable z. To obtain a Green function g°''\y^z) for 
A(s) — satisfying (4.21) above we seek a function 

fiR{a)3z^g''\y,I) 
which is a solution to the Helmholtz equation, 

(A(s)-A«)^«''(-,?) = 0, 

defined for z G Dij(a) with z G Y>fi[a)\{y} (recall that pr(i) = z and pr(y) = j/), 
with monodromy multiplier e"^'^*^^ in a counterclockwise circuit of the branch point 
aj. Also, to obtain a Green function we want g°''''{y,z) to differ from g^^siy^z) by 
a smooth function when y is close to z. Note that, because of the transformation 
properties we want for g°''''{y, z) it will have singularities on the set 

V := { pr-i(a;) x : x G Djt(a)} 

as well as on the diagonal of Di?(a) x Di^(a). 

To make this precise we let (p^z) denote a Cq° function on with 

(j){z) = 1 for 1^1 < 1 
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and 

(piz) = for \z\ > 2, 

and define 

Then for y fixed and different from a j for j = 1 , . . . , n and e small enough the 
support of z ^ (p^ (z — y) will not contain any of the branch points {ai, 02, ... , Cbn}- 
In such circumstances it is easy to see that that the function, 

2; ^ 0e 

has support in an elementary neighborhood of the point z = y. An elementary 
neighborhood, U, of a point y G Dij(a) is one for which the inverse image of U 
under the projection, 

pr : Di^(a) T)R{a), 
splits into components with the property that 

is a diffeomorphism. The sets Vq, are indexed by homotopy classes, a, of paths that 
join the base point ao to y. Using this property we will now indicate how to obtain 
a lift of the function, 

U xU 3 {x,y) ^ (j)e{y - x)gi^^s{x, y), 
to a function defined on 

PT~\U) X pi-\U)\D 

with prescribed multipliers. Let ao denote the base point in Dij(a) as above. Let 
ao denote a path joining ao to y representing fixed homotopy class [ckq] . Once [ao] is 
fixed we can index the components of pr~^(C/) by homotopy classes of loops based 
at y. Furthermore the choice of allows us to identify the homotopy group of loops 
based at ao G T)R{a) with the homotopy group of loops based at y G Di?(a). Thus 
if we let n denote the representation of the fundamental group of (Dij(a), ao) which 
results from assigning the number e"^'^*'^ to the generator which starts at ao and 
makes a simple counterclockwise circuit of aj not winding around any of the other 
points ttk for k 7^ j, then this naturally gives a representation of the fundamental 
group of (Dij(a), J/) as well; namely, 

7r[7] := Trfag ^7Q!o] for [7] G 7ri(Dij(a), y). 
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Now we define a fift of 

9l,six,y) ■= My - x)gK,s{^^y) 

which is defined on U x U\{{x,x) : x e D_R(a)} by 

9l,si^,y) ■■= ^]g^^,s{^,y) for x e V[^,«o] and y e V[^2ao]- 

The function y) is defined on pr~-'^(t/) x pr~^(C/)\r', and [jk] G 7ri(Dij(a), y) 

for /c = 1, 2. It is not hard to see that this extension is independent of the choice of 
[cto], for if ctQ is another path from gq to y one has 

and 

7r[Q;(^"'^Q;o]7r[Q;o~^Q;o] = 1- 

It is now natural to extend g^^giv^ ^ Djj(a) x Djj(a)\D by setting it equal 

to on the complement of pT~^{U) x pT~^{U)\'D. Now fix y G Djf(a). We will 
construct g°''''{y, T) so that up to a smooth error it looks like gf. ^(y, 2) whenever z, 
the projection of z, is close to y, the projection of y. More precisely, we require that 

{X{s)-A,){g^^\y,I)-glMI)} (4.24) 

be free of singularities as a function of ? e Dii{a). When z is not in the fiber over 
y = pr(y) we want (A(s) — A^)g"''\y,^ = 0. Thus with the formal cancellation 
of the delta functions in (4.24) the only remaining contribution comes when the 
differential operator part of A« hits the function (f)^{z — y) in g^^g{y,T) (note that 
the multiplicative part of which does not contain any derivatives is absorbed in 
the delta function that arises when hits g^^^siy-, z)). Now choose e small enough 
so that 2e will produce an elementary neighborhood U above. Then (f)2e{z — y) is 
identically 1 on the support of (p^{z — y). Using this one can easily see that there 
is an explicit second order differential operator, p{dz, d^), whose lowest order terms 
are first order in dz and dz and for which we want 

(A(.) - A,) {g^\y,T,-glM^} = {p{dz,dz)M^-y)}gl%{y,T,. (4.25) 

Since p{dz,dz) kills (^^(2; — y) near z — y it follows that the right hand side of 
(4.25) is a smooth function of 'z in C^^. We will now show that for any / e 
whose support projects onto an elementary neighborhood it is possible to find a 
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multivalued function F, transforming as / does under deck transformations, so 
that 

(A(s)-A«)F = /. (4.26) 

Once we have this existence result we can then recover g"''''{y, i) from the solution 
to (4.25) above by adding in g^ g{y,T). 

Now we will recall the Hilbert space technique of §3 for finding a solution to 

(4.26) . Let Ha,i denote the Hilbert space completion of with respect to the 
norm derived from the inner product 

{F, G)= f \KkF{I)KkG{I) + m2F(l)G(F)| dfx{z). 

Note that the integrand descends to a function on Y)ji because the monodromy 
multipliers of F and G cancel. Now suppose that / G and the support of 
/ projects onto an elementary neighborhood. For v e consider the linear 
functional 

I mv{T, df,{z) :=< f,v>. (4.27) 

Since (/i, h) > < h,h >, for all h e C^i, it is clear that the linear functional 

(4.27) is continuous on the Hilbert space Ha^i- Thus by the Riesz representation 
theorem there exists an element F e Ha,i so that 

<f,v>= {F,v). 

Since the formal adjoint of is one sees from this last equation that F is a 
distribution solution to (4.26). Recall from the discussion in §3 that F will be locally 
in the Sobolev space (at least away from the branch points) and so by standard 
regularity results for elliptic equations, one can conclude that F must actually be 
C°° away from the branch points aj for j = 1, 2, . . . , n. 
We now use this result to solve the equation 

ix{s)-A,)F^{p{d,,d,)Mz-y)}9l:s{y,^- (4.28) 

Let g"-'^{y,z) denote the sum of the solution F{y,z) to (4.28) and gl.^g{y,z). Then 
the reader can check that g""'^ satisfies the following three conditions: 

Gl: The map z — > g'^'^y, i) is a C°°(Dij(a)) solution of the Helmholtz equation, 

{X{s) - A^)g-'\y,I) = 0, 
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loi z ^ y and z ^ aj for j = 1, 2, . . . , n which transforms under the deck transfor- 
mations Rj as follows, 

G2: If t/ is an elementary neighborhood of the point y then in the complement 
Di?(a)\pr~^(f/) the functions, 

z^g^^y,^, 

z^Kk9'''\y,I), 

z^K*,_,g'^'\y,I), 

have absolute squares that are integrable with respect to the natural hyperbolic 
measure, diJ,{z), on D^. 

G3: For some elementary neighborhood, U, of the point y the difference 
is smooth as a function ofzE pr~^(t/). 

Next we want to recall the consequences of G2 for the local expansions of ? — > 
g'^'^y, i) about the points aj for j = 1, 2, . . . , n. Because g"''''{y, i) is a multivalued 
solution to the Helmholtz equation with monodromy 6"^"^*^^ about the point aj it 
follows that fir"''' (y, T) has a local expansion, 

nez 

valid for z near aj. Recall now that Ij is chosen so that < Ij < 1 and that the 
leading singularity in vi^K{z,a) and a) goes like r' in both cases, where r is 

the distance from z to a. It is now a straightforward exercise to verify that the 
condition G2 implies the following restriction on the local expansion for g"-'\y,'z), 

oo oo 

f\y. ^ = 5] ai{y)vn-i,A'^^ a,) + PiivX+i.A'^^ «^)' (4-29) 

n=l n=0 

for z near aj. We will see that this restriction on the local expansions of g"''\y, z) 
together with Gl and G3 and the fact that z — > fif"''(y, 5) is square integrable on 
IDr suffice to characterize g"''''{y,z). We also wish to characterize g°''''{y,z) as a 
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function of the first variable y. Since A„ is Hermitian symmetric one might expect 

that 

g-^\y,T> = g-'i{I,y). (4.30) 

Observe that the right hand side of (4.29) transforms under y — > Rjy by the mul- 
tipher e^'^*'^ which is appropriate for g^-^^y^z). Also since gK,s{yiz) = gK,,siz,y) it 
follows that the right hand side of (4.30) differs from the Green function g^^siUiZ) 
along the diagonal by a smooth function — this is additional evidence that the right 
hand side of (4.30) should agree with g°'''-{y,T). We will now confirm (4.30) by 
a Stokes' theorem calculation along the lines of the standard argument for such 
symmetry. Now choose x and y two distinct points in Dij(a) and let 

g{Tj = g'''\l,y) 

and 

fiT, = g-'\l,x). 

Substitute these last two functions in (4.7) and integrate the resulting equality over 
the complement of the union of the balls of radius e about each of the points aj and 
the balls of radius e about both x and y. The left hand side gives since 

{X{s) - A,)/ = (A(s) - A,)^ = 

on this set and Stokes' theorem reduces the integral on the right to an integral over 
the circles of radius e about the points a-, for j = 1, . . . , n and x and y. Using the 
results for the local expansion of g'^'''{z, x) given in (4.28) above one finds that each 
of the integrals, 

/ 'fi{f,Ki^g)-ip{K^f,g), 

vanishes in the limit e — > 0. Finally, using the fact that g°''\x,y) differs from the 
Green function for A(s) — A^ by a smooth function near the diagonal one finds that 
the contribution to the integral of 

vQ^K^g) - (p{KJ,g) 

over the circles of radius e about x and y is just 

9{x) - f{y) =g°"\x,y) -g'^'^iy^x) 

in the limit e — > 0. This last quantity must vanish and this establishes the Hermitian 
symmetry of the Green function we are considering. 
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The Green function for the Dirac operator in the presence of branch 
points. Next we turn to the construction of a Green function for the Dirac operator 
acting on by using formulas (4.18) and (4.19) with Qi^^g replaced by g"''''. This will 
not give us the Green function we want because the local expansions will not have 
the "minimal singularity" type that we desire. However, by subtracting suitable 
wave functions from the results we are able to pass to the desired Green function 
in much the same fashion that we constructed the wave functions Wj from the 
canonical wave functions Wj. In analogy with (4.18) and (4.19) we now define 
g'^'i by 



Gi'\y,z) = (m - Dk) 



S2\y,z) = (m - Dk) 













(4.31) 
(4.32) 



where we have explicitly noted the dependence of (7"^!, on («;, s) to avoid confusion. 

We now record the local expansion results for (4.31) and (4.32) which are obtained 
from (4.29). There are two different cases to consider depending on whether < 
< I or i < < 1. In case < < ^ we write Xj — Ij and the local expansions 
for Q"'''' are about z = aj are 



n>- 



(4.33) 



n>i 



where all the sums are over elements of Z + i and we have written «);^fe(J, a; m) and 
wIj^(z, a; m) for (2.22) and (2.23) in order to recognize explicitly the dependence on 
the parameters k and m (it is —m which occurs in (4.33)). 

In case ^ < < 1 we write Aj = Ij — 1 and the local expansion results about 

z 



aj are 



Gi\y, ^ = X] ^n,i^n-\^,k{z, aj- -m) + J2 (^i,iK+Xj,k{^, aj; -m), 

n>-i 



(4.34) 



n>-i 



n>| 



In (4.33) and (4.34) we have suppressed the dependence on parameters y, /c, m, and 
Xj for the coefficients ct^ ^ and ^ to avoid further complicating expressions that 
are already notationally overburdened. 
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To obtain the Green function we desire our object will be to eliminate those 
terms in the local expansions (4.33) and (4.34) with n < |. We can do this by- 
subtracting suitable multiples of the wave functions Wj{—m) and Wj{—m) without 
introducing additional singularities on the "diagonal" y = z. As is the case in 
(4.33) and (4.34) it is now convenient to explicitly identify the m— dependence of 
the wave functions Wj and W* since the Green functions Qj' {y, i) are solutions to 
the equation {m-\- Dk)il^ = in the F variable instead of the original Dirac equation 
(m — Dkji^ = 0. As always our notation for the response functions will suppress 
dependence on parameters that do not play a role in the current calculations and 
will rely on characteristic names for the parameters to make our intentions clear to 
the reader. Thus Wj{z, —m) is our short hand notation for Wj{z, k, A, — m). Define 

G^'^ {y, ^ = {y,^- Yl ^ -™) • 

' * 2 ' 

j:\j<0 

As the reader can easily verify the subtractions in (4.37) remove all terms in the 
local expansions on the right hand side with n < | and we find 

Gi'^{y,z) = ^^e^^/y)wn-xA^^ aj-, ~m) + f^/y)w*^+^.{z, ttj-, -m)^ , (4.38) 

where we have once again omitted the k dependence of wi^k and Wij^ for brevity. 

We turn next to a characterization of G'^'^{y, z) along the lines of the charac- 
terization Gl, G2, and G3 of g''^\y,z). 

Proposition 4.0 The components G"'^(y, ?) of the Green function constructed 
above for (m — D^) acting on are uniquely characterized by the following three 
properties: 

(1) For i = 1,2 the vector valued functions z — > G'^''^{y,z) are solutions to the 
equation, 

{m + Dk)G't'\y,T> = 0, 

when z is not equal to aj for j = 1, . . . ,n and not equal to y. Under the deck 
transformation Rj, G"' (y, T) transforms as follows, 

G^'^(y,i?,i)=e-2-^^G^'\y,i). 

(2) About each of the points aj for j — 1, 2, . . . , n the functions G"'^(y, i) have 
restricted local expansions given by (4.38) above. In the complement of a compact 
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(4.37) 



neighborhood of the set {y, ai, 02, . . . , an} the function z Gl'^iy^ z) has a norm 
which is square integrable. 

(3 ) For some elementary neighborhood U, of y the difference, 

Gt'\y,^-Gtiy,T,, 
is smooth as a function ofze pr~^{U) and G\{y, z) is defined by, 



Gliy,T, 

Gliy,^ 



Proof. These three properties of G^' (y, z) are simple translations of the properties 
Gl, G2, and G3 for the Green function g"''\y, z). The fact that these three condi- 
tions characterize G"' (y, z) can be understood as follows. Suppose that G"' (y, z) 
and i) both satisfy all the conditions in Proposition 4.0 and consider the 

difference 

Aiy,^ = G^'\y,^-H:^\y,^. 

Then (3) of Proposition 4.0 implies that A{y,z) is free of singularities near z — y 
and (1) and (2) imply that z ^ A(y, i) is a solution to (m + Dk)A{y,'z) = with 
monodromy multipliers e"^'^*^^ about aj and which is globally in L'^{dii). However, 
the restricted local expansion (4.38) implies that for each j = 1,2, ...,n either 
Cj{A) = or c*(A) = and this together with Theorem 3.0 implies that A = 0. 
QED 

The derivative of the Green function G^'^. Next we will identify the low 

order expansion coefficients e\ .{y) and f{ .{y) in (4.38) as the components of wave 

2 2 '* 

functions for the Dirac equation. Suppose that f{z) is a solution to the Dirac 
equation (m — Dk)J = which transforms under the deck transformations Rj as 
follows, 

/(i?,i)=e-2-^V(i), 
and which has restricted local expansions, 

/(^ = {^n'^ri-Xj (z, aj; m) + KK+Xj (z, a^; m)} . (4.39) 

n>-i 
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To lighten the notation for the next calculation we write 

Now let D^{y) denote the disk of radius e about y, and let C^{y) denote the circle 
of radius e about y. Then for z in the complement of the union, 

D,{y)UD,{ai)U---UD,{an), (4.40) 

we have 

Dkf{I)=mf{I), 

Dug{z) = -mg{z). (4.41) 

If we now make use of (4.41) in (3.5) we find that 

d \ —^-^^r-pridz idz > = 

1 _ Ml 1 _ N_ I 

in the complement of the set (4.40). If we now integrate this last equality in the 
complement of the set (4.40) and use Stokes' theorem to rewrite the result as a 
circuit integral on the boundary, then we find 

/ I '^^^^19 idz ^^^^idz \ — — f < —^-^^^idz idz \ . 

UJo.ia^Al-^-^ 1-^ J 4(.)ll-¥ J 

(4.42) 

Of course, this makes sense since f^gj is a well defined function on Dij(a)\{y}. We 
now wish to evaluate both sides of (4.42) in the limit e 0. As in Theorem 3.0 it 
is useful at this point to introduce the local coordinate Uj = T[—aj]z in the C^[aj) 
integral and the local coordinate u = T[—y\z in the Ce(j/) integral. Then in geodesic 
polar coordinates (r^, 6j) and (r, ^) we have 

u 

and the asymptotics are given for small Vj and small r. Making use of this asymp- 
totic parametrization in (4.42) and replacing some factors arising from the change 



Re'hh 



Rr 
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of variables by their limits as e — > 0, one finds 

lim ^ei? / {7i^2K)e-''^ + /s^iKOe^'^ } dO^ 

' (4.43) 
= - lim eR r {lrg2{u)e-'' + J^g^{u)e''] dO. 

We have taken the liberty of scaling the result by a factor of 2. Next we use (4.38), 
the asymptotics of the wave functions wi and given in (3.1) and (3.2), and the 
local asymptotics of the Green function for z near y to evaluate both sides of (4.43) 
in the limit e — > 0. We find 

E ^-^^ {^ei ,(^ +?7/i ,(^| = -im. (4.44) 

^ ' m H L 2 2'* 2 2'* J 

7 = 1 



where we used 



r(A)r(i-A) =^^"^"^^- 

We now use (4.44) to evaluate the low order expansion coefficients e\ .(y) and 

2 '* 

/i . (y) in terms of the wave functions Wj and . First choose 

2 

then 



a 



=0, 

2 

and (4.44) specializes to 



fl,{y) = -'^W.{y,m)i, (4.45) 



where Wi/(- • ■)i is the z*'* component of Wi^{- ■ •). Next choose 

f{I) = W:{z,m), 

then 



ail =0, 

2 

2 
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and (4.44) specializes to 



el/y) = -^W*{y, m)i. (4.46) 

Now we are ready for the principal result of this section. We write G'^j^{y^z) for 
the j*^ component of the vector G"'^(y, 5). We are interested in calculating the 
derivative of G"j^(^, z) with respect to the a variables. As in §2 it is convenient to 
introduce the scaled variables 

, at, 
^^ = ^- 

Taking the derivative of G"'^(y, F) in the parameter bi, kills the singularity on the 
diagonal y — z. Thus the function 

I ^ d,^G^^\y,T, 

will be a solution to (m + Dk)ip = with monodromy multiplier e~'^'^^^^ associated 
with the deck transformation Rj. This function will be determined uniquely by its 
local expansions at aj for j = 1,2, ...,n. Using (2.32) and (2.34) to differentiate 
the local expansions (4.38) we find that 

dKGt'\y,T, = -(1 - |&.|2)-ie| .(y)W-.(i',-m). (4.48) 
Using (2.33) and (2.35) to differentiate the local expansions (4.38) we find that 

dt^G^'\y,T> = -(1 - K\Y'fliiy)W:i^, -m). (4.49) 
Finally substituting (4.45) and (4.46) in (4.48) and (4.49) we find 



db.G'^-\y,^ = ^W^j(y,m),iy.(i',-m),-, 



(4.50) 
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§5 Deformation Equations 

A holonomic system. In this section we will find deformation equations for the 
low order expansion coefficients of the response functions, 

= Sj^W_i+)^. + J2 {<uf^n+Xj + KuK-x^ } (5-1) 

n>0 

and 

W* = Sj^w*_i_^, + |4^«;n+A,. + di^wl_^. } . (5.2) 

n>0 

It might be helpful for the reader to observe a parallel between the developments 
here and one way of understanding the Schlesinger equations from the theory of 
monodromy preserving deformations of linear equations in the complex plane. In 
the Schlesinger theory the fundamental object is a multivalued analytic function 
z a) which takes values in the general linear group. The matrix valued 

function Y{z, a) has fixed (i.e., independent of the aj) monodromy matrices Mj in 
simple circuits of the branch points aj and satisfies an ordinary differential equation, 

dz ^-^ z — aj 

j=i -I 

in the z variable. This equation is the analogue of (5.3) below. The fixed mon- 
odromy condition allows one to extend this differential equation to a differential 
equation, 

' ^(--''.) ' 

in the a variables. This is the analogue of (5.20) below. The compatibility conditions 
for the z and a equations gives the Schlesinger equations, 

daAy = - ^-^^^^^^d{aj, - ap), 

which is the analogue of (5.26) below. We derive the deformation equations (5.26) 
in a manner precisely analogous to this. One difference is that the internal com- 
patibility conditions for the z and a equations do not give anything new in the 
Schlesinger case. In our case there is some interesting information that needs to be 
extracted from the internal consistency of (5.3). 
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Our first tool will be the following linear relations, that arise when the infinites- 
imal symmetries Mj are applied to Wj^ and W* : 

biMiw, + M2W, = Yl {^-i^^i^ + /^-M^;} . 

In these relations the matrices e, /, h, a, P, 7, and S are matrices that depend on 
the branch points aj for j = 1, . . . , n but not on z. In the case of e, /, g, and /i we 
have also partly anticipated the form that the relations will take. This will simplify 
some later results and will avoid the necessity of introducing further definitions. The 
terms —\Wy and \W* are stuck out in front of the first two equations instead of 
being combined with e^^ and h^^^ to avoid introducing the Kronecker 5^^. The Jj^^ 
that appears in the last equation is not the Kronecker delta — the possible confusion 
with this will shortly be addressed when a, /3, 7, and 5 are all eliminated in favor of 
e, /, g, and h. 

These relations (5.3) will be deduced using a combination of (2.20) and (2.32) 
— the results that describe the action of the infinitesimal symmetries Mj on local 
expansions and the results for centering the operators Mj at different points in the 
disk. As the reader can easily check using (2.20) and (2.32) the linear combinations 
on the left hand side of (5.3) are all chosen so that the leading singularities in 
the local expansions cancel at level n = — |. Since VF^ and W* for ji — 1, . . . ,n 
are a basis for the wave functions whose local expansions contain terms no lower 
than n = — I and which are in LF' near |2;| = i? relations of the form (5.3) follow 
immediately. To find the coefficient matrices e, /, and etc. one need only compare 
the local expansions on both sides of (5.3) at level n — One finds, 

e= (A; + A) + [ai,5] , 
/ = (5bi5-bi), 
^ = (ci - BciB) , 
h={k + \)- [di,S], 
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(5.3) 



(5.4) 



where we used the notation 





A* 

= 

n- 


1 

2 


,.(1- 








n— 


1 

2 


.(1- 


IV 






= c„_ 


1 
2 


.(1- 


l&M 






n- 


1 
2 


.(1- 


IV 





(5.5) 



Also we have written B for the diagonal matrix with uv*^ entry given by bi, and 
{k + A) for the diagonal matrix with entry {k + X,^) on the diagonal. The analogue 
of the display (5.4) for the coefficients a, etc. is 



[ai,i?' 



a = ~2B{k + X- i) 
p = {h, - B%^B^) , 
7 = (ci - B^CiB^) , 
5 = 2B{k + X+l)-[duB\ 



(5.6) 



From (5.4) and (5.6) one may deduce some simple relations among the matrix 
coefficients, 

a = B - eB - Be, 
p=-f-BfB, 
-i = g + BgB, 
6 = B + hB + Bh. 



(5.7) 



Next we will deduce some less obvious relations from the compatibility requirements 
of (5.3). Consider the Lie algebra relation, 



(M3M2 - M2M3 - M2)W = 0, 
where we've written W for the column vector, 

Wi 



(5.8) 



W = 



We now treat (5.3) as equations for M2W and M^W in terms of MiW and W and 
W* and as equations for MiW* and M^W* in terms of M2W* and W and W*. If 
we now use (5.3) to eliminate the appearance of M3W, M2W, M3W*, and MiW* 
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in favor of MiW, M2W*, W and W* one finds that the coefficient of MiW in the 
resulting equation is, 

[e,B^] + [a,B], 

which vanishes as a consequence of (5.7). The coefficient of M2W* in (5.8) is, 

{B^fB^-f + BpB-p)B, 

which also vanishes as a consequence of (5.7). Since W and W* are linearly inde- 
pendent we can equate the remaining coefficients of W and W* to in (5.8) and 
one finds, 

a+[e,a] + 2Be + Bp^ + B^fB^-PgB-B = 0, 
ep + 2BfB + Bp6 + B^fB6-afB- ph = 0. 



In a similar fashion if one starts with the Lie algebra relation, 

(M3M1 - M1M3 + Mi)W* = 0, (5.10) 

and uses (5.3) to eliminate MiW* , M^W*, M2W, and MsW in favor of MaM^*, 
MiW, W and W* then one finds the coefficient of MiW is, 

(B''gB^-g + ^-B^B)B, 
which vanishes as a consequence of (5.7) and the coefficient of M2W* is, 

[h,B'] + [B,6], 

which also vanishes as a consequence of (5.7). Equating the coefficients of W and 
W* to in what remains of (5.10) one finds, 

/17 + SgB + 76 + B'ja - 2BgB - B^gBa = , 

_ _ — —2 — (^-11) 

S - [h, S] + ^fB + BjP - 2Bh - B gBp -B = 0. 

Both (5.9) and (5.11) simplify if one uses (5.7) to eliminate a, (3, 7, and 6 from these 
equations. Substituting (5.7) in (5.9) one finds 

[B,e^-fg] = (5.12) 

and 

X - BXB = 
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where 

X = ef- fh. 

However, X — BXB = 0, implies X = so the second consequence of (5.9) can be 
written 

ef-fh = 0. (5.13) 

Observe also that (5.12) implies that the matrix — fg is diagonal. In a precisely- 
analogous way one can substitute (5.7) in (5.11) to get 

[h^-gf,B] = (5.14) 

and 

X - BXB = 

where 

X = ge- hg, 

and as before this implies 

ge-hg = 0, (5.15) 

and (5.14) implies that the matrix /i^ —gf is diagonal. The equations we've found in 
this fashion do not determine the diagonal parts for — fg and h'^ —gf - However, we 
can determine these by looking at the coefficients of w'i+Aj local expansions 

of the first and third equations in (5.3). One finds for the first equation, 

1 _ \h.\2 "^2l-2 + ^jJOji' H 1 _ |, 12 + + 2)^lu 



l-|6,f "i^ ■ l-\b, 



If we set v = j this simplifies to. 



2' 



/i=l /i=l 



-6jm2(-| + Aj) + (/e + A^- + 1)0-1^. = ^ ej^a|^ + ^ /?>&mc(m' ^^'-^^^ 

where we have taken the liberty of collecting some of the terms involving a\ . on 
the left hand side. If we look at the coefficient of wij^^. in the third equation of 



(5.3) then set 1^ = 3 and collect some of the terms involving a\ . on one side then 

2 J 

we find, 

(1 + |6,f )m2(A,- - i) - 2h^{k + A,- + l)a\. 

^ ■ ^ - ■ (5.17) 



It 

/i=i /i=i 
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If one now multiplies (5.16) by 2bj and adds the result to (5.17) then one finds 
(after a little simplification), 

(e^ - fg),, = e + ^. 
Since we already know that — fg is diagonal it follows that 

e^-f9= (k^ + I- (5.18) 

In a precisely analogous fashion one can use the second and fourth equations in 
(5.3) to show that 

h'-gf-ik' + '^^i. (5.19) 



The Holonomic Extension. Next we consider extending the differential equa- 
tions (5.3) to the a variables. To do this it is convenient to use (2.32)-(2.35) to 
compute the local expansions for the exterior derivatives dhWu and d^W*. One 
finds for the local expansion of d^Wi, about the point aj, 

(1 - ) dbW^ ^^-i+Xj [-^j'^dbj] 



+ W_i 



Sj^k_{Xj,bj)-a{^dbj 



I * 

+ w 1 . 



-b\ db, 



+ wi^^. da\^ + a{^k+{Xj, bj) - aijbj - 5j„m2{\j - \)dbj 



I ^ 

+ 



+ bik_{Xj,bj) - h'^dbj 



+ 



where we've introduced the abreviation 



k±{\j,bj) = {k + Xj ± ^){bjdbj - bjdbj). 
For the local expansion of d^W* about the point aj one finds, 

{i-\bj\')d,w:=w*_s_. [s.^db,] 



+ w_ 



2 -^3 



5j^k+{Xj,bj) - d{dbj 



T+Xi 



+ Wi . 



+ Wl 



+Xi 



—c\ dbn 

d{d\^) + d\^k-{Xj,bj) - di^dbj - 5j,,mi{-Xj - i) 
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Using these results one can easily check that the terms at level n — cancel in 
the local expansions of 

and 

dbw: +dKM2wt. 

Thus one has relations of the form, 



(5.20) 



By comparing terms in the local expansions of both sides of (5.20) at level n = \ 
one finds that, 



1 — 

F = -dBhiB^ -hidB , 



(5.21) 



G = - dBciB^ - cidB, 
1 — \B\ 

The deformation equations. The combination of the Dirac equation, the first 
two equations of (5.3) and (5.21) determine an extended holonomic system for 
the response functions W and W*. We will deduce deformation equations for the 
coefficients e, /, g, and h by examining the compatibility conditions between (5.20) 
and the first two equations in (5.3). The two compatibility conditions we will 
examine are 

{Msdb - dbMs) W = 0, 
{M^db - dbMs) W* = 0. 

We use the equations in (5.3) and (5.20) to express all the terms in the preceeding 
equations in terms of MjW, MjW*, W and W* for j = 1, 2 and = 1, 2, . . . , n. In 
a single term in each of the equations this also requires the use the commutation 
relations for the Lie algebra of Mi, M2, and M3 and one should also make the 
observation that the most singular terms containing M^W and M|VF* manifestly 
cancel. The condition that the lowest coefficients in the local expansions of the 
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resulting expressions must vanish leads to the following identities: 



[E, B] - [dB, e] = , 

c/S/# + fdB + BFB - F = , 

[B, H] - [dB, h]=0, 



(5.22) 



dBgB^ + gdB - BGB + G = . 



All these identities are simple consequences of (5.21) and (5.4) but it will be useful 
to record them here for future use. Once one has identified these relations all the 
remaining terms involving MjW and MjW* can be written as linear combinations of 
W and W* using (5.3). Since W and W* are linearly independent we can equate the 
coefficients of the resulting equations to 0. We find the first form of the deformation 
equations, 



We will first eliminate the appearance of ct, /3, 7, and 6 from these equations using 
(5.7) and then we use the identities (5.22) to rework the equations as explicit equa- 
tions for df and dg. We will illustrate the calculations involved for e and / and 
leave the analogous reductions for g and h to the reader. 
First substitute ^ = g + BgB into (5.23) to get 

de = -{dBfB + BF)g - {dB0 + BFB)gB + [E, e] + FgB - fBG. 

Now use (5.22) to substitute F - fdB for dBfB^ + BFB on the right hand side of 
the last equation. After some cancellation one finds 



difB) 
d{gB) 



dh 



de 



-{dBfB + SF)7 + [E, e] + FgB - fBG , 
- {dBfB + BF)5 + {EfB - eF) + {Fh - fBH) + F , 
(BG - dBgB)a + {Ge - gBE) + {HgB - hG) - G , 
(BG - dBgB)P + [H, h] + {GfB - gBF) . 



(5.23) 



de = f{dBgB - BG) - {dBfB + BF)g + [E, e]. 



It will prove useful to introduce, 



F = - 



{dBfB + BE) = dBhiB + Bh^dB 



(5.24) 



and 



G = dBgB -BG = dBciB + BcidB. 



(5.25) 
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The equation for de becomes 

de = fG + Fg + [E,e]. 

Now we work on the equation for df. Again we start by replacing S on the right 
hand side hy B + hB + Bh. One finds that, 

{df)B + fdB = F(B + hB + Bh) + EfB - eF + Fh - fBH + F. 

Next we use (5.22) to replace FB by fdB — F in this last equation and make some 
cancellations to find, 

{df)B = FhB + fdBh + EfB - eF - fBH + fhdB. 

Using (5.22) we can replace the term f{dBh — BH) that appears on the right of 
this equation by 

f{hdB - HB) = efdB - fHB, 
where we used fh = ef. One finds, 

{df)B = FhB + EfB - fHB + e{fdB - F). 

But using (5.22) we can replace fdB — F in this last result by FB and then cancel 
the common factor of B on the right to obtain 

df = Fh + eF + Ef- fH. 

Doing precisely analogous calculations for g and h one finds, 

de ^fG + Fg+[E,e], 

df =eF + Fh + Ef - fH, ^ ^ 

^ ^ ■' •' ' (5.26) 
dg =hG + Ge + Hg-gE, 

dh =gF + Gf +[H, h]. 



Further identities. There are a number of identities among the low order expan- 
sion coefficients of the response functions W and W* that are associated with the 

transformations m < m and (k, X) <— {—k, —A) which we will now record. First 

observe that, 

wi,k{-,m) = wi^k{-,-m) , 

(5.27) 

f^lki-^'m) = wlj,{-,-m). 



1 



1 
-1 
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Because we are interested in the dependence on the parameter m we have made 
it exphcit in the preceeding equation when we did not always do so before. From 
these relations it follows by comparing local expansions at level n = —\ that, 



-1 
1 

1 
-1 



Wi,{k, A, m) = Wi,{k^ A, —m) , 
W*{k,\m) = W*{k,X, -m). 



Now comparing local expansions in (5.28) at level n — \ one finds, 

a^iik, A, m) = at (k, A, -m) , 
bt(k, A, m) = -bt(k, A, -m) , 

2 2'^ 

ct(k, A, m) = -ct(k, A, -m) , 

2 2 

c?^ (/c. A, m) = dt (k, A, -m) . 



(5.28) 



(5.29) 



Since the dependence of the low order coefficients on the parameter m is so simple 
we will surpress the explicit notational dependence on m with the convention that 

for X = a,h,c and d. If bj^{k, A, —m) (or c^j,{k, A, — m)) occurs we will replace it 
with —bj^{k,X) (or —c^^{k,\)) without further comment. 

For k — Q the Dirac operator is a real operator. One consequence of this is 
that there is a natural conjugation that acts on the A; = solutions to the Dirac 
equation. This conjugation is. 



w 



1 

1 



w. 



(5.30) 



It is not difficult to check that this conjugation acts on the local wave functions 
(2.22) and (2.23) as follows. 



1 

1 

1 

1 



w/,fc(-,a) = w*i a) , 



(5.31) 
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Using (5.31) to compare local expansions at level n — —h one finds, 



1 

1 

l" 

1 



w,{k,x)^w:{-k,-x), 



(5.32) 



W*{k,X)^W,{-k,-X). 



Comparing local expansions at level n = | in (5.32) one finds, 



d'i(k,X) = al(-k,-X), 



cUk,X)^bU-k,-X). 



(5.33) 



Finally we will note some identities that follow from the derivative formula for the 
Green function (4.50). To make use of this formula it is convenient to modify the 
Green function with (5.32) in mind. We introduce 



8 



m 



2i? 



1 

1 



and (4.50) translates into 



1 



s.(l-|6.P) 
1 



W^{-k, -A, m) (g) W^{k, A, -m) , 
W*{-k, -A, m) (8) W*(k, A, -m) 



(5.34) 



If we now use this equation and the local expansion formulas for dbW,^ and d^W* 
given above to compute the local expansion coefficient of 



w_i_^.{-,ai)(S)w_i+^.{-,aj) 



in the equation. 



one finds after some simplification 

s~^ai {—k, —A) — a 1 {k, Xys~^ — diagonal. 



(5.35) 



where ai is the matrix with i/fx matrix element a^^ and X'^ denotes the transpose 
of the matrix X and s denotes the diagonal matrix with i/v element Si, = sin(7rAi/). 
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In a similar fashion if one calculates the local expansion coefficient of 

in the identity, 
then one finds 

s'^bi {-k, -A) - ci {k, Xys-^ = 0. (5.36) 

Now the combination of (5.35) and (5.33) produces the following relation between 
the deformation variables h and e above (in which the parameters are understood 
to be {k, A, m)), 

h = s{l- \B\'')e*s-\l - \Bf)-\ (5.37) 

where e* is the conjugate transpose of e. 

The combination of (5.36) and (5.33) produces the relation, 

(^s-^6i(/c,A))* =s-^6i(/c,A), (5.38) 

from which it follows that 

{fs{l-\B\')y = fs{l-\B\'). (5.39) 

One also has 

(^s-^ci(A;,A))* = s-^ci(A;, A), (5.40) 

from which it follows that 

{s-\l - iBlY'gy = s-\l - \B\^)-'g. (5.41) 
We summarize these developments in the following theorem: 

Theorem 5.0. The low order expansion coefRcients e, f, g, and h deGned by (5.4) 
above satisfy the non-hnear deformation equations, 

de =fG + Fg+[E,e], 
df=eF + Fh + Ef-fH, 
dg =hG + Ge + Hg-gE, 
dh =gF + Gf+ [H, h] . 
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The coefEcients e, /, g, and h also satisfy the relations, 



ef 


-fh = 


0, 


ge 


-hg = 


0, 




-f9 = 






-9f = 





where — k'^ + ^ . If we further introduce, 

f = /s(l-|S|2), 
g = s-\l-\B\Y'9, 
then we have the additional symmetry relations, 

f* = f , 

* 

h = sil-\B\^)e*s-\l-\B\'')-\ 



We have one final observation to make concerning these deformation equations. If 
we combine the commutator equation ef — fh = with the final expression for h 
we find, 

ef - fe*. 

If for some reason / turns out to be invertible (this will happen in an interesting 
special case we consider later on) then we can use the non-linear relations above to 
eliminate g and h from these equations in favor of e and /. In particular one finds, 

h = r'ef 

g = r\e^-M\ 

Theorem 5.1. If / is an invertible matrix and we define, 

F = Fs(l-|Sp), 

^ (5.42) 

G = s-\l-\B\^r'G, 

then the deformation equations for e and f are 

de = fG + Fg+[E,e], 

^ ^ (5.43) 
df = eF + Fe* + Ef + fE*. 
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The matrices e and f satisfy the conditions, 



ef = fe*, 

f * = f , ^ ' ' 

with g determined by e and f , 

g = f-i(e2 - M^) = e*f-^e - M^f-i, (5.45) 

and G determined by g. 

We leave it to the reader to confirm that (5.43)-(5.45) incorporate all the informa- 
tion about e, /, g, and h that we have found so far in the event / is invertible. 

The two point deformation theory in a special case. In this subsection we 
will integrate the deformation equations in a special case of the two point problem. 
Our principal result is that it is possible to integrate the deformation equations in 
terms of a Painleve transcendent of type VI. 

It is useful to start by determining what the deformation equations have to say 
about the action of the rotational symmetry bj — > e'^^bj. This corresponds to the 
infinitesimal symmetry, 

d d ja-d\ (h - d 

rot = © • • • © -^rr- = I bi— bi^ © • • • ©^ bnTr, bn 



ddi d9n \ dbi dbj \ dbn "db, 

One sees immediately that 

rfS(rot) = iB , 
dB{rot) = -iB. 
It is then a simple matter to calculate, 

E(rot) = -ie + i{k + A) , 
F(rot)=z/5, 
G(rot) = -igB, 
H{Tot) = -ih + i{k + X). 

The deformation equations become 



n 



(ie(rot) = i[k + A, e] i[X, e] 
df (rot) =i[X J]. 
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From these equations it follows that, 

eie'%u . . . , e^'b^) = e'^'e{bu . . . , b^)e-'^' , 
f{e'%,, . . . , e^%n) = e'^'fib,, . . . , 6„)e-^^^ 



(5.46) 



Next we want to consider what happens for the infinitesimal version of the SU(1,1) 
one parameter group, 

cht sht 
sht cht 

The vector field associated with the action of this one parameter group on the disk 
of radius R is 

V = Vi ® ■ ■ ■ ®Vn 



with 



It follows that 



1 

R 



d_ 
db, 



1 

'r 



_d_ 



1 



dB{v) = -{l-B'), 



dB 



One computes. 



E{v) 
F{v) 



1 
1 



-[{k + \){B-B)--{B + B)+eB + Be] , 



R 



if + BfB) , 



1 



G{v) = --{g + BgB) , 



H{v) = 



1 

R 



and 



{k + X){B -B)- ^{B + B)-hB- Bhj , 



1 



F{v) = --{fB + Bf) , 

G{v) -^i9B + Bg) . 
A little computation using the relations among e, /, g, and h shows that 

de{v)^^[X~^(B + B),e] , 

df{v) = ^[X-l(B + B)J]. 



(5.47) 
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Since the matrices 



A 



are diagonal they commute for different values of the matrices B and it follows that 
the differential equations for e and / may be explicitly integrated along the orbits 
of the one parameter group action, 



t 



cht 
sht 



sht 
cht 



for p e Djj. We wish to parametrize the two points 61 and 62 in the following 
fashion, 



hi 



b2 



ri 



^,,^ch_(^25 + sh_^ 
^ sh(^)5 + sh(^) 



where 



The parameters for bi and 62 are 0, ^i, r, and ri. Note that (r, ^) are essentially the 
geodesic polar coordinates for 62 relative to 61. The reason for this parametrization 
of bi and 62 is as follows. First rotate both hi and 62 by —9i , then take the result 
of this and flow along the appropriate v— orbit by — ^. The result is 



61^0, 62 



S. 



Finally rotate and 6 simultaneously by —9. One finds that 

bi ^0, 62 ^ th (0 . 

It follows that the symmetry fiows that we have determined above suffice to deter- 
mine e(6i, 62) and 62) for all values of (61, 62) in terms of e(0, t) and /(O, t) for t 
real. Next we will use the deformation equations to write down ordinary differential 
equations for e (0, t) and / (0, t) . Let 



d_ 
dt' 



Now we substitute (0, t) in the deformation equations and evaluate them along the 
vector field u. One finds, 

^0 0" 
t 



B 



B 










t 



(5.48) 
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and 



dB{u) = 
dB{u) = 





1 



1 



(5.49) 



To make use of the symmetry conditions to reduce the number of variables it is 
desirable to work with e and f rather than e and /. Indeed at this point we will 
confine our attention to the case in which Xj > for j = 1,2. In this case the 
response functions Wi, are in and the matrix —s~^bi{k,X) is the matrix of 
the inner products of these reponse functions and is consequently positive definite. 
Using the fact that \bj\ < 1 one easily sees that this implies f is positive definite 
(and hence invertible). Hence the deformation equations become 



de = fG + Fg+ [E, e] , 

di = eF + Fe* + Ei + fE% 



with 



g = e*f-^e-M2f-\ 
and G determined by g. At this point it is useful to introduce 



A; + Ai e+ 
e_ k + X2 



and 



Kch ij) esh -0 



(5.50) 



(5.51) 



(5.52) 



(5.53) 



where c > 0, k, and if) are real and ee = 1. This parametrization incorporates the 
known diagonal for e and the fact that f is a symmmetric positive definite matrix 
which is related in a particular fashion to the matrix of the inner product of 
wave functions. The last algebraic relation ef = fe* implies, 

(ee_ — ee_)shi/j = 0, 
(ee+ — ee+)shi/; = 0, 
{k~^{1 — t^)e+ — K,e-)ch.'ip — Xeship = 0, 

where A = A2 — Ai. Now define 
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If we suppose that sh'^ ^ then we find that ee+ and ee_ are real and 

^+ " 2(1 -t^) *^^^^^^^^' 
e_ = ^(^-AthV^), 



where C — Now write 



A/ ji — /t ~|~ A j , 

;U = /ci + /c2 = 2/c + Ai + A2 , 
A = /c2 - ^1 = A2 - Ai , 



for J = 1, 2 and 
and 

Then one easily calculates, 
Eiu) 

F{u) 

Giu) 



f = 



fii in' 

^21 f22 

gll gl2 

g21 g22 



2t(A;2 



1-^2 
fi2 





'0 


0" 




' 


e+ 





1 


e_ 






1372 122 

-gl2 



F{u) = 
G{u) = 



T^&22 _ 

2t 



2t 





f22 



g22 



If we now substitute these into the equation for df{u) we find 

df{u) = N 



where 



iVn = -t-i(e+f2i + e+fi2), 
t + t-^ 



N12 
N21 
N22 



t + t-^ 



e+f22 - t ^e_fii + 



2t(A;2-|). 



12 



22 



-t-i(e_fi2 + e_f2i) 



l-t2 

t e_iii H -5 — 121 



2t 



1-^2 



f22- 



(5.54) 
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Next we convert this using (5.53) and (5.54) to find, 



iVl2 

N22 



= —CK, 



—ce 



(^ + AthV')shV': 



r 



—ce 



1-^2 



1-^2 



■ship 



-CK"4"^(^ - AthV')sh'0 - CK~^2tchi;. 



Write f = % and the equation df{u) = N for c' , e' , k' and ip' becomes, 

c~^cch.il) + K~^Kch.'i(; + ship ip' = Nn , 

c~^csh.ip + t~^es\vip + chip ip' = c~^e~^Ni2 , 
c~^c'shijj — e~^e'sh'0 + chip ip' — c~^eN2i , 
-2tchV' + (1 - t^) (c-^c'chV' - K~^K'chip + ship ip') = c~^nN22 ■ 

The very last equation simphfies somewhat when combined with the expression for 
N22- One finds. 



c ^c'chip — K ^ k' ch ip + sh ip ip' 



r 



1-t- 



:(^ - AthV')sh'0. 



Solving these equations for c ^c', k ^k', e ^e' and ip' one finds. 



1-^2 



0, 



Ar 



iP' = 



l-t2 



(a^ - i)t 

l-t2 



ch'0 sh'?/;. 



Next we work out the differential equation associated with e. Substituting the 
results for F and G into the differential equation for e above, one finds. 







2t 



72&22 







1-t 

]-z|2-fl2g22 
— lZ:|2-f22g21 



2f f 
1372122 



S+[E,e] 
+ [E,e]. 



(5.55) 
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We also calculate 



[E,e] 



A + (/i-l)t2 



Combining (5.55) and (5.56) we find, 

2t 



-e+ 
e_ 



e' = 



2t 



fl2g22 - ^(i_^2) e+ ' 

A+(/X-l)t^ 
f22g21 H — e_ 



Use 



to obtain, 



1 - t2 t(l - ^2) 

g = f-i(e2-M^) 

(1 - t2cli^?/')cg2i = -e (A;i - + e_e+) shV' + «;//e_chV' 
(1 — t^cli^'0)cg22 = — e/ie+sli'0 + k(/c2 — + e_e+)cli'0 . 



(5.56) 



(5.57) 



(5.58) 



Now substitute this last expression for g2i in the second equation in (5.57) along 
with the parametrization for f22 and the expression for e_. Then multiply both 
sides of the resulting equation by k, make use of the equation for k~^k' given above 
and the observation that 

e'-A'th^ 

e_e-|- ^ 



4(1 -t2) 



to obtain. 



lei.. ^ , ,x A(l -th^V') + (a*- 1)^^ ^ , ,x 

tchV^shV' /,9 ,o „,,o f^-A^th^V * .r.^ 
t/ich tp 



1 - t^ch^V 



As a check on this equation we remark that the variable '0 in the Poincarc disk has 
the same significance as the variable il) in [8,12,15]. Introducing t = ^ and letting 
— > oo one finds that the limiting form of (5.59) is 

(r^\ =A2(l-th2V;)thV'+^sh2V' 
or \ or 2 
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which is equation (3.3.46) in [12]. In the special case A = this equation was first 
found in [8]. 

Next consider the substitutions 



One finds, 



2 1 

S = t , W = T — 

chV 



di/j 1 di/j 1 dw 

ds 2t dt 2w\J\ — w ds 

1 ^s{l-s)^ + s{^i-l){l-w)). 



\/l — w V ds 



Now convert (5.59) to a differential equation in s by multiplying both sides by 
and substituting s for t^. In the resulting equation make the substitution 
w = l/ch^ip and multiply both sides by w\^l — w. One finds, 

wVl-w4- 7^ ) (s{l -s)^ + {{fM- l)s - Xw){l - w) 

ds Wy/l — w J \ ds 

+ ^ (^(1 - + ((/^ - 1)^ - - ^) 

Xw + (u — l)s ( , .dw ,, X ^ X/ X 1 .X 
^ ' s{l-s)— + {{ii-l)s-Xw){l-w)) (5.60) 



2s(l - s) V ds 



tu-s V 2(1 -s) 

^ '^5(l-s)^ + s(/x-l)(l-tt;) 



w — s \ ds 
It is simple to compute, 

- d 1 2 — 3w dw 



wy/l — w 



ds w^/l — w 2w{w — 1) ds 



Substituting this last expression in (5.60) and simplifying the resulting equation 
one finds, 

2 \w w — 1 w — s J \s s — 1 w — s J 

w{w-l){w-s) f{l-4M'^)s(s-l) {n-lfs_^X^^ ^^'^^^ 



s2(l-s)2 y 2{w-s)^ 2w^ 2 
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This is Painleve VI (see, e.g., [4]) 



with 



1 (\ 



w 



+ 



1 



+ 



1 



2Vm' w — 1 w — s 



+ 



s — 1 w — s 



w 



w{w — l){w — s) — 1) 7(s — 1) /3s 



a = 



P 



2 ' 



7 = 0, 
, 1 - 4M2 



We summarize these developments in the foUowing theorem: 

Theorem 5.2 Suppose that Xj > for j = 1,2. Let bj = ^ and consider e(&i, 62) 
and 62) as functions of the scaled variables bj. Define A = A2 — Ai and /j, = 
2/c + Ai + A2. Write 

^A; + Ai e+{t) " 
e-{t) k + X2 



and 



e(0,t) = 
f(0,t) = c(t) 



K(t)chV'(t) e(t)sh?/'(t) 
e(t)shV'(t) -t^jchV'lt) 

Tien tie algebraic relation ef = fe* implies, 

en 



-2(1 -t^) 



(^ + AthV^) 



e_= — (e-AthV'), 

for a reai valued function ^ and e of absolute value 1. The deformation equation for 
f implies, 



(//- l)t ,2 , 



l-t2 



= 0, 



Ar 



t2 



thV, 



t 



-1 



l-t2^ ' l-t2 
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chi/; sh'^. 



The last equation allows us to eliminate ^ in favor of •0 and if:' . Once this is done 
and the further substitutions, 



s = t^ , 
1 

w 



arc made, the deformation equation for e becomes the type VI Painleve equation 
(5.61 ) above for w. 

§6 The Tau Function 

The Grassmannian formalism. In this section we will introduce the r— function 
for the Dirac operator with branched singularities whose Green function was de- 
scribed in §4. Suppose that S is an open subset of T)r with a boundary dS that 
consists of the union of smooth simple curves in D^; (we assume that the closure 
of S is contained in the open set Dr). Suppose that / is a smooth function in the 
Sobolev space H^CDji) which satisfies the Dirac equation 

(m - Dk)f{x) = for X e S. 

Then since the integral operator associated with the Green function inverts the 
Dirac operator we have 

fji^) = ("^ - Dk)f ■ Gj{x,y) diJL{y) for x e S. 

Jr>n\s 

The integral can be confined to HrXS since (m — Dk)f{y) = for y e 5. Since the 
Green function Gj{x, y) is a solution to 

{m + Dk)yGj{x,y) = 

in the second variable when x ^ y, it follows that the expression for fj{x) can be 
written 

fj{x)= I Urn - Dk)f{y) ■ Gji^) - f{y) ■ (m + Dk)Gj{x,y)\ 

which becomes, using (3.5) and Stokes' theorem, the following boundary integral 
representation for fj{x), 



/,(.)=2/ ^^^i(^#M,^_^i^(^^,,,W (6.1) 



9S\ 1-^ 1-^. 
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in which each component of the boundary, dS, is given the standard counterclock- 
wise orientation. If / G H^{dS) is prescribed arbitrarily then we can also regard 
(7.1) as a formula for the projection of / onto the space of boundary values on dS 
of H^{S) solutions to the Dirac equation in the interior of 5". The complementary 
subspace for this projection is the space of boundary values on dS of solutions to 
the Dirac equation in H^{'Dfi\S). 

Our principal tool in the discussion of the tau function for the Dirac operator is 
a formula analogous to (6.1) with the Green function G replaced by G"''^. To explain 
the significance of this formula it will be useful to describe the transfer formalism 
[12] and its relation to "localization" away from singularities. It is natural when 
considering the Dirac operator with branch type singularities at the points aj for j = 
1, 2, . . . , n to localize the operator away from the branch cuts. One natural way to 
do this in the hyperbolic disk is to transform the disk into the upper half plane by a 
fractional linear transformation and then to draw horizontal strips around horizontal 
branch cuts emanating from each of the points aj for j = 1, 2, . . . , n. There are many 
ways to transform the disk into the upper half plane and it is clear that, provided the 
branch points aj are all distinct, one can choose a suitable such transform so that 
the second coordinates of the transformed branch points are distinct. Supposing 
this to be done we can then choose the strips Sj containing aj to be pairwise 
disjoint by making them sufficiently narrow. It is certainly possible to transform 
the Dirac operator into upper half plane coordinates and work exclusively in the 
upper half plane, so that the description of the localization is geometrically simple. 
However, because we wish to use the many formulas which we've written down in 
the disk but have not written down in the upper half plane this is not economical. 
Instead we adopt the expedient of using the geometrically natural terminology in 
describing the localization via horizontal strips in the upper half plane but when 
doing calculations we transform the strips Sj back into the disk where they become 
crescent shaped objects with the sharp ends meeting at a common point on the 
boundary of the disk. The reader should not find it difficult to keep in mind that 
any reference to horizontal branch cuts or strips implicitly assumes that one has 
chosen an appropriate upper half plane coordinate system. 

Until now it has been convenient to work on the simply connected covering 
Dj^(a) of Dj^(a). However, at this point it is simpler to draw horizontal branch cuts 
(rays), £j, emanating to the right from each of the branch points aj and to work 
with functions on T>ii{a)\{ii, i^, . . . , in} that have appropriate branching behavior. 
We will now describe the localization of the singular Dirac operator, D"^'^, which 
we wish to consider. Let S = USj denote the union of the strips Sj. In the exterior 
of the set 5" the operator, D"'^, will act precisely as the ordinary Dirac operator 
acts. However, in order for this operator to adequately mirror the singular Dirac 
operator whose domain contains functions that are branched at the points aj with 
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specified monodromy we restrict its domain. We will now define the subspaces of 
H^{dS) that are important in our description of localization of the singular Dirac 
operator D"^'^. We will say that a function F is locally in the null space of D"^'^ at 
ttj provided that for some e > the function F has a local expansion, 

^(^) = {an{F)wn+Xj {x, aj) + bn{F)wl_^.{x, a^)} , 

valid for x e Bf:{aj), the ball of radius e about aj. In this description it is understood 
that a branch cut ij has been chosen for aj and that branches for F , Wn+x- , and 
'^n-Xj fixed. Suppose now that Sj is a horizontal strip containing aj in its 

interior. We define a subspace of H^{dSj) in the following fashion, 

Definition. A function f on dSj will be in the subspace Wint{aj) provided that it 
is the boundary value of a function F deGned in Sj with the following properties: 

(1) F is a branched solution to the Dirac equation in Sj with branch cut £j and 
monodromy e^'^*^^ . 

(2) F is locally in the null space of D"''^ at aj. 

(3) F e H^{Sj\{B^{aj)uej}) for some e > 0. 

If each branch point aj is contained in a strip Sj and the strips Sj are pairwise 
disjoint then we also write, 

Wintia) = I^int(ai) e I^i„t(a2) • • • W^int(an), 

to define a subspace of H^{dS). Next we define a subspace of H^{dSj) comple- 
mentary to W\^x.{aj). 

Definition. We define the subspace VFext(%) to consist of functions f defined on 
dSj which are the boundary values of functions F defined in 'Df;\Sj which satisfy 
the Dirac equation (m — Dk)F = outside Sj and such that F G H^{'Dii\Sj). 

Note: For both of the preceeding definitions the functions of interest are in the 
first Sobolev space near the boundary and so the Sobolev embedding theorem 
implies that the boundary values are assumed in the norm on the boundary. 

Finally we let Wext denote the subspace of boundary values on dS of functions F 
defined in Dr\S which satisfy the Dirac equation (m — Dk)F = in the exterior of 
S and such that F e H^{'Dji\S). Note that this subspace is certainly not the direct 
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sum of the subspaces VFext(oj). The foUowing formula defines a projection that is 
fundamental for us, 




Theorem 6.0 The subspaces Wintia) and W^^t cLre transverse to one another in 
H2 (dS). The operator P"''^ defined in (6.2) above is the projection on VFint(o) along 

The proof of this theorem is parallel to the calculation used to establish (4.7) in 
reference [12] and so we will not repeat it here. 

There is another projection that will be important for us. Let Pj denote the 
projection of H^{dSj) on Wint(aj) along VFext(oi)- We define 

F{a) = Pi e P2 ® • • • ® Pn- (6.3) 

We are now prepared to discuss the determinant bundle formalism that will allow 
us to define a r-function for the singular Dirac operator D"-'^. Fix a point a° = 
(a5, 02, . . . , a°) with 0° a1 if j 7^ k. Fix a collection of pairwise disjoint open 
strips Sj with e Sj. Choose e > so that for each j = 1, 2, . . . ,n the ball of 
radius e about a° is contained in Sj. That is, 

P,(a;)c5,. 

Let H = H^{dS) and write 

= W^int(a°)©Wext (6.4) 

for a distinguished splitting of the Hilbert space H. Write Pq for the projection of H 
on Wint{aP) along W^^t- We are interested in the Grassmannian, Gro, of subspaces 
of H which are close to VFint(o°) in the following sense. A closed subspace W is in 
Gro provided that the map, 

Po:W^^W^int(a°), (6.5) 
is Fredholm with index 0, and the map, 

(/-Po):W^^Wext, (6.6) 
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is compact. Except that we have additionaUy specified the index of the first pro- 
jection this is essentially the definition of the Grassmannian that Segal and Wilson 
consider in [16] and the reader can find a detailed theory of such Grassmannians in 
the book [14] by Segal and Pressley. The restriction to index in (6.5) means that 
we confine our attention to the connected component of the Grassmannian, Gro, 
containing Wint{a^)- It is simpler to discuss the det* bundle over the connected 
component of the Grassmannian and it will suffice for our purposes. The first result 
we require is 

Theorem 6.1 If aj e B^{a^) then Wint(a) G Gro. 

The proof of this theorem is precisely parallel to the proof of Theorem 4.1 in [12] 
and so we refer the reader to this paper for details. 

The localization of the singular Dirac operator we are interested in can now be 
described as follows. The operator, D"''*^, acts as the ordinary Dirac operator in the 
exterior of the union of strips S and its domain is the subspace of H'^{Gr\S) whose 
boundary values belong to the subspace VFint(«)- We define a determinant for this 
family of Dirac operators by trivializing the det* bundle over the family of subspaces 
Wirxt{o) £ Gro. This trivialization is then compared with the canonical section for 
the det* bundle and the result defines the r— function for the Dirac operator. 

Recall that an invertible linear map F : VFint(a°) — >^ is an admissable frame for 
the subspace W G Gro provided that Po-^ '■ W^int(a°) Wintia^) is a trace class 
perturbation of the identity. We will now introduce two difi'erent frames for the 
subspace VFint(o). The first such frame will define the canonical section of the det* 
bundle and the second will be used to provide the trivialization of det* over the 
family or subspaces 

a Wint{a) e Gro. 

First we show that the restriction P"--^ : VFint(o°) W^int(o) inverts the projection 
-Po : Wint{a) VFint(o°)- Prom this it follows that the restriction of P"''-^ to VFint(a°) 
is an admissable frame for VFint(a) which defines the canonical section of the det* 
bundle. The argument is simple. Suppose that w G Wi^tia'^) and w = Wa + 
with Wa G VFint(a) and We G Wext- Then P'^'^w = Wa and writing Wa = w — We 
we see that Po^a = w. Evidently, this is just an expression of the fact that the 
complementary subspace Wgxt is the same for both projections Pq and P°^'^. 

The second frame arises from considering the direct sum of "one point" pro- 
jections: 

Proposition 6.2 The restriction of F{a) to the subspace Wint{a^) is an admissable 
frame for Wint(a)- 
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Once again this is identical to Proposition 4.2 in [12] and we refer the reader to that 
paper for more details. 

We now recall that the fiber in the det* bundle over a subspace W G Gro can 
be identified with equivalence classes of pairs (w^a) where w : Wint(«'^) W is 
an admissable frame and a is a complex number. The equivalence relation which 
defines the fiber is {wi,ai) = {w2, if and only if 

ai — a2det{w2^wi). 

In this representation the canonical section of the det* bundle is given by 

Gvo3W ^ {w, det{Pow)) e det* 

where w is any admissable frame for W. Since PqP"''^ is the identity on VFint(a°) it 
follows that we may regard 

Gro9 W^i„t(a)^(P"'^|o,l) 

as a representation of the canonical section a{Wint{a))- We use the notation P'*''^|o 
to signify the restriction of the projection P"^'^ to VFint(fl°)- We now use F{a) to 
define a trivialization of the det* bundle over the family of subspaces VFint(o). Define 

(^(W^int(a)) = (F(a)|o,l). 
We may then define a determinant T(a, a°) for the Dirac operator as follows, 

= ^ deto(P"'^lo 'F(a)) = deto(F(a)|o-^P"'^)-^ (6.7) 

In this formula deto refers to the determinant of an operator restricted to the sub- 
space Wint{a^)- The second form of the determinant as a reciprocal is included 
since it will be slightly simpler for us to calculate the logarithmic derivative of the 
T-function in this form. We can now state one of the principal results of this paper: 

Theorem 6.3 The r— function defined by (6.7) has logarithmic derivative given 

by, 

dlogr = J2 i_\b.\2 + : (6-8) 

where the coefRcients a\. and d\. are the local expansion coefRcients for Wj and 
W* found in (3.22) and (3.25) above. 
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Proof. For simplicity we write P{a) = P"'-^ in the following calculation. We also 
write d for the exterior derivative with respect to {ai, 02, . . . , a„}. The rule for 
differentiating determinants when applied to the final term in (6.7) above yields, 

dlogT(a,a°) = -Tvo{d{F{a)-^P{a))P{a)-^F{a)), 

where Tro is the trace on VFint(a°) and it is understood that both F{a) and -P(a) 
are regarded as maps from VFint(o°) to VFint(fl) in this last formula. However, we 
know that P(a°) inverts P{a) restricted to Wint{0'^)- Thus 

P(a)-^F(a) = P(a°)F(a) 

when both sides are restricted to Wint{a^)- In a precisely similar fashion F(a°) 
inverts F{a) restricted to Wiat{a^)- Thus we find, 

dlogr(a,a°) = -Tro(F(a°)d(P(a))P(a°)F(a)). 

But P{a){I — P{a'^) = so that d{P{a))P{a^) = dP{a) where dP{a) is now regarded 
as a map on all of H. Combining this observation with the fact that the range of 
F{aP) is all of Wintia^) we can remove the subspace restriction on the trace to find, 

dlogT(a,a°) = -T:(F(a")(dP(a))F(a)) = -Tr(((iP(a))F(a)F(a")). 

Since F{a){I - F{a^)) = we find that F{a)F{a^) = F{a). Thus 

dlogT(a, a°) = -Tr{{dP{a))F(a)). (6.9) 

The reader might note that although r(a, a*^) depends on its logarithmic deriva- 
tive does not! Since it is clear that T(a°, a°) = 1 from the definition, the fact that 
the logarithmic derivative is independent of a° implies that there exists a function 
T(a) such that 



The function r(a) is the function which we would like to identify as a correlation 
function in a quantum field theory. This will have to await further developments. 
However, one should observe that 

dlogT{a, a°) = dlogT{a). 

We will now use the derivative formula (4.50) and the formula (6.2) for the projec- 
tion to evaluate the trace (6.9) in terms of the low order expansion coefficients for 
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the wave functions Wj and W*. The derivative formula (4.50) and the projection 
formula (6.2) together imply, 



TDt \ m^RW^ix.m) f W*{y,-m)^ 



-12) / — \:^My)^dy, 

v\ ) Jds 1 - 



(6.10) 



and 



^^^^ ^ 'rn^RW*{x,m) f W^{y,-m)^ 



-h{y)i'dy 



4^.(1 -|6.P) 1-%^ 



(6.11) 



The appropriate trace is, 



-TV(F(a)a..P(a)) = - ^^^^^^^ X, ^ZW ^ ^ 



i?2 



m'^R f W*{y,-m)^F{a)Wt,{y,m)i . ^ 

z dy, 



^sAl-MJas 1-^ 

(6.12) 

and 



4.^(1 - |6.P) ^5 



4..(l-|6.|^)ya5 1-^ 

(6.13) 

According to (3.5) both of the integrands in (6.12) and (6.13) are exact away from 
the singularities at aj for j = 1, 2, . . . , n. The contours dSj can be closed about 
the singularities aj and the same asymptotic calculation that leads to (4.4) above 
shows that, 

-Tr{Fdt^P) = ] .^ a{FW,)l (6.14) 

and 

-Tr{Fd,^P) = —L_b{FW:)l, (6.15) 
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where we used the notation, 

n 

for the local expansion coefficients of g about the point a,^. Note that the contours 
dSj with j ^ u do not make a contribution in this calculation. To compute FWi, 
and FW* on dS^, we need only subtract from W^, and W* the "one point versions" 
of these two functions (just the translate of the appropriate multiple of the function 
to '^j^)- Doing this we find, 

aiFW,)", =aK^ 

2 2 

and 

b{Fw:n = di^. 

2 2 

This finishes the proof of Theorem 6.3. QED 

We will now provide the connection between the coefficients in the formula for 
the logarithmic derivative of r given in (6.8) and the deformation theory of §5. If 
one collects all the terms in (5.16) that involve a\ . on the left hand side using 

2-? 

Cjj = k + Xj then one finds, 

a\. = bjm2{Xj - I) + XI ^iM^p + X -^iA^/^' ^^""^^^ 
The analogous calculation for d\ . yield, 

d{. = bjmi{-Xj - I) - X ~ J^djf^^f^^^^- (6-17) 

The right hand sides of both (6.16) and (6.17) are completely determined by the 
deformation variables e, /, g and h. 

The log derivative of r for the two point function. In this subsection we will 
write out the formula for the logarithmic derivative of the tau function in terms of 
the Painleve transcendent of the sixth kind that was found in the integration of a 
special case of the deformation equations described at the end of §5. We suppose 
therefore that Xj > for j = 1,2 and that 6i = and 62 = t. In this special case 
(6.8) becomes, 

|logr(0,t) = ^{a|,, + 44. 
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From (6.16) and (6.17) one finds, 

a| 2 = tm2{X2 - I) + ^^loi i + */22c| 2 

and 

d 1 2 = tmi{-X2 - |) - ^2lC^| 1 - ^^22^1 2- 

Expressing the off diagonal elements of a, b, c and d in terms of the deformation 
parameters one finds, 

a^i , =t-^(l -t^)ei2, 

^1,2 = - /22, 
4,2 =^22, 

4 1 = -t-i(i-t2)/ii2. 



Thus, 



a| 2 = tm2(A2- ^)+t ^(l-t^)e2iei2 + t/225'22 



4,2 = tmi{-X2 - I) + t-'(l - t^)h2lh^2 + t/22^22. 

Now we observe that 

m2(A2 - I) = mi(-A2 - \) = +k^-{k + A2)' 

and since 

g = .-i(l-|S|2)-i^, 
f = /.(l-|S|2), 

it follows that 

f 22921 = f22g22- 

Eliminating h in favor of e* one has, 

h,2="-^{l-tY'el,, 
h2i = — (1 -^^) 621 , 

Si 

so that, 

h2ihi2 = 621612 = 621612, 



where we have used the fact that 621612 = 6_e+ is real. Thus 
|logr(0,t) = ^ 



4 logT(0, t) = {2tm2(A2 - I) + 2t-\l - t')e_e+ + 2tf22g22} • (6.18) 



89 



Using (5.53) (5.57) and (5.58) one finds, 



-efie+ship + K{k2 — + e_e+)chi/; 



g22 = 



= K- 



1 - t2ch> 

-fi^^sh^ + (kl -M'+ '^i^) ch^ 
1 - t^chV '' 



so that, 



r22g22 



-2/;(4 + Atliv)riiv shv + ((2 - A2tlrv)cli%- + 1(1 - f-)(A:| - i\/2)ch> 

4(1 - t^ch^V) 



As a check on (6.18) we can, as we did for (5.59), introduce t — ^ and then take 
the hmit — > oo. One finds that the hmit of cilogT(0, -^) as i? — > oo is 

1 f^^y _ ^2^^2^ _ ^2^2gi^2^ I 

This should be compared with (4.5.42) in [15] and (5.39) in [12]. The overall sign 
difference with (4.5.42) is a consequence of the fact that (4.5.42) is a formula for 
the log derivative of the Bosonic tau function in SMJ. 

Next we make the substitutions, w = l/ch^-i/' and s = in the expression for 
f22g22- Substitute the result in (6.18) along with (5.54) for e_|_ and e_. Eliminate ^ 
in the result with the substitutions 

CthV' =^ (s(l - s)w' + s(l -ii)iw- 1)) , 
w'^{l — w) 

Noting that 

we find after some calculation 



a, s{s-l) fdw w-lV 

— logr 0,v^ = . , ^ r 

OS Aw[w — l)[w — s) \ds s — 1/ 



w — s 

fi^ A2 ^ 4M2 - _ 



(6.19) 



A{s-l)w 4s(s-l) 4s 4(1 -s) 
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The formula on the right hand side of this last equation bears some striking resem- 
blance to the formula in Okamoto [7] for a Hamiltonian associated to the Painleve 
VI equation. It does not, however, appear to be a Hamiltonian for Pyi for any- 
simple choice of canonical coordinates. 

It would be interesting to study the behavior of this tau function as t ^ 0, 1. 
The asymptotics for t near should match onto the asymptotics of the analogous 
tau function in the Euclidean domain. The behavior as t — > 1 would follow if 
one knew the appropriate connection formula for Pyi- We hope to return to this 
question elsewhere. 
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